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Abstract 

The master equation describing non-equilibrium one-dimensional problems 
like diffusion limited reactions or critical dynamics of classical spin systems 
can be written as a Schrodinger equation in which the wave function is the 
probability distribution and the Hamiltonian is that of a quantum chain with 
nearest neighbor interactions. Since many one-dimensional quantum chains are 
integrable, this opens a new field of applications. At the same time physical 
intuition and probabilistic methods bring new insight into the understanding of 
the properties of quantum chains. A simple example is the asymmetric diffusion 
of several species of particles which leads naturally to Hecke algebras and q- 
deformed quantum groups. Many other examples are given. Several relevant 
technical aspects like critical exponents, correlation functions and finite-size 
scaling are also discussed in detail. 
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1 Introduction 



Our understanding of nonequilibrium statistical physics is far behind that for the 
equilibrium theory. Even simple models may pose a formidable problem if one wants 
to approach them analytically. In this paper, we shall consider two different types of 
such problems: the diffusion-limited chemical reactions and the critical dynamics of 
classical spin systems. It will be shown that in one dimension these problems can be 
mapped onto quantum chain problems which are often integrable and on which a lot 
of progress was recently achieved [0. As a result, many new predictions concerning 
the nonequilibrium statistical physics problems follow. 

The study of the diffusion-limited chemical reactions has stimulated a vast amount 
of research since the first investigation of Smoluchowski many years ago [[J . Examples 

are given by the bimolecular reactions, A + B # C + D where the two species A and 

B diffuse and react to form the two new species C and D and k and g are the forward 
and backward reaction rates, respectively. The simple case of irreversible reactions 
for which g = 0, C is a inert product and D is not present has been extensively 
investigated since the original work of Zeldowich [H. Despite their simplicity those 
systems exhibit a very rich dynamical behaviour. For homogeneous initial conditions 
and in low dimensions, the diffusion mechanism is not efficient enough to mix the 
particles. As a result a spatial segregation occurs and accordingly, a reduced number 
of reactions between the two species is possible. This results in a slowing down of 
the dynamical evolution called anomalous kinetics M. The evolution of the system 
is not properly described by the usual rate equations, since the fluctuations play a 
crucial role. Another interesting situation is when the reactants are initially separated 
in space; then a react ion- diffusion front is formed during the evolution ||. Here 
again the properties of the front are drastically influenced by the fluctuations in low 
dimensions ||. We shall restrict ourselves to the homegeneous case here. Several 
approaches have been used to study such systems: numerical simulations |7J], scaling 
arguments || and analytic arguments based on the theory of the Brownian motion ||. 
The main results obtained concern the decay of the number of particles of one species 
U and the temporal evolution of the gap developing between the particles of the two 
species due to the segregation process flQ| . However, the analytic results are scarce. 



Very little is known concerning the behavior of quantities like the space and/or time 
dependent two-particle correlation functions. Accordingly, new theoretical approches 
allowing the computation of such quantities are desirable. 

Let us first introduce several models which have been studied in the literature: 



1. The coagulation model: WL\ one considers molecules of one species (say A), that 



diffuse in a milieu and react as: 

A + A^A (1.1) 



2. The annihilation model: jl2| the A molecules diffuse and annihilate 

A + A^dl (1.2) 
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where denotes an inert state which decouples completely from the dynamics. 

3. The two-species trapping reaction: two types of molecules A and B diffuse 
and one of them, B, is "trapped" by A: 

A + B^A (1.3) 

4. The two-species annihilation reaction: || two types of molecules A and B diffuse 
and annihilate 

A + B (1.4) 

We are interested in the long-time behaviour of systems like the four examples 
given above. Considering quantities like the mean concentration of A particles c^(t), 
in general we expect the following two types of behaviour, as t — > oo 

c A (t)~{ r \, (1.5) 
I exp — t/T 

where a is some constant and r is known as relaxation time. Throughout this paper, 
we shall refer to the first type as "massless" or "critical" , while the second case will be 
denoted as "massive" . This terminology is borrowed from field theory and equilibrium 
statistical mechanics. All four reactions defined above have a critical behaviour in 
the sense that for a given initial concentration of particles, the long time behaviour 
of the concentrations has an algebraic fall-off. For example, in the reactions (|1.1| ) and 
( |1.2| ) and in one-dimension, the concentration c^(t) of particles behaves like 

c A {t)~t- 1 ' 2 (1.6) 



One can allow for reversible (or back) reactions, corresponding to g ^ 0. In the 
long time limit (t > a local equilibrium state is reached |14|]. One is then in a 



"massive" regime, in which the relaxation towards equilibrium is exponential. 

A different type of nonequilibrium problem is the one of critical dynamics of 
classical spin systems. Let us consider for example the case of a classical Ising model 
on a lattice. The system is prepared in an initial (nonequilibrium) state. One would 
like to know how fast the equilibrium state (Gibbs state) will be reached. When not at 
(static) criticality, the system relaxes exponentially towards equilibrium. In general, 
the relaxation time r scales with the spatial correlation length £ as the temperature 
approaches its critical value 

r~f (1.7) 

where z is the dynamical critical exponent. However, counterexamples will be given in 
this paper where £ diverges but r remains finite (no critical slowing- down!) As these 
classical spin systems do not have an intrinsic dynamics, the dynamics is thought 
to come from the interactions between the spins and a heat bath modelling the fast 
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degrees of freedom not included in the classical Hamiltonian. The dynamics of these 
models is thus formulated in terms of a master equation for the probability that a 



spin configuration is realized at time t 15 1. Several cases have to be distinguished 



depending on the presence of macroscopic conserved local quantities. The simplest 
case is the purely relaxational one in which there are no local conservation laws and we 
shall restrict to this case throughout the paper. For example, the initial magnetization 
relaxes towards its equilibrium value. All the physics is put into the transition rates 
appearing in the master equation. Several choices are possible compatible with the 
condition of detailed balance, which is the condition insuring that the stationary state 
will be the Gibbs one. Despite its apparent simplicity, the solution of this problem 
is very difficult and little can be said analytically even in one dimension |l6fl . The 
mapping onto a quantum chain Hamiltonian, which will be explored throughout the 
paper, will turn out to be a very useful tool to clarify some controversies present in 
the literature pf, fl8|. 

Both types of nonequilibrium problems (diffusion-limited chemical reactions and 
critical dynamics) can be described in terms of a master equation for P({f3},t), the 
probability that a configuration {/?} of the systems is realised at time t. It turns 
out that it is suitable to map this master equation problem onto a quantum chain 
problem |16 , [19], 2(J. The corresponding equation of motion reads 



d t P({(3},t) = -HP (1.8) 

where H is directly related to the transition rates appearing in the master equation. 
This will be detailed in Section 0. The question of knowing what kind of dynamical 
behaviour has the model (power law or exponential relaxation) amounts to know in 
which phase (massive or massless) of the phase diagram we are. 

Moreover, in reaction-diffusion processes H turns out to be non-hermitian and 
has often the following particular structure 

H = H + H 1 (1.9) 

where Hq is a known integrable Hamiltonian with a larger symmetry than H (for 
example, in two-states models it can be the XXZ quantum chain in a Z field). Hi 
is non-hermitian and has a lower symmetry but does not affect the spectrum of H . 
Thus, if Hq is massless, it follows that H is massless as well. In problems of critical 
dynamics, H has again often the structure of Eq. ( |1.9| ) but with a new meaning. H 
is now hermitian, Hq is again a known Hamiltonian (it can be the same as the one 
occuring in reaction-diffusion processes) but Hi is now a perturbation term, if we 
approach criticality (small temperatures). 

When writing this paper we were faced with two problems. The first one was 
that we realized that we are left with many more open questions than answers. This 
is kind of nice because we hope that this is an invitation for other people to look 
closer at the subject. The second problem concerns the pedagogical presentation of 
the paper since it addresses two different communities: people doing nonequilibrium 
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statistical mechanics and who are familiar with the physical problems treated in 
this paper but not with integrable systems and those familiar with the Bethe ansatz 
and two-dimensional field theory but not familiar neither with the physical problems 
discussed here nor with methods of computing nonequilibrium averages which are 
different from those techniques used in equilibrium statistical mechanics (vacuum 
expectation values). We thus suggest two approaches to this paper. One for the 
"mathematician" , the other for the "physicist" . 

The "mathematician" should start with Appendix C (the last chapter of the pa- 
per). There we remind the reader of the definition of the Hecke algebra which depends 
on a parameter f3 = q + q~ l (the significance of q will become apparent immediately). 
As is well known, if a quantum chain can be written as a sum of generators of a Hecke 
algebra, through Baxterization [^1 one can associate to the chain an integrable ver- 
tex model. In this paper, we will consider only chains with 2 L and 3 L states ( L is the 
length of the quantum chain). Accordingly, we are going to look for some quotients of 
the Hecke algebra. To various quotients of the Hecke algebra one can associate a rep- 
resentation given by the (m/n) Perk-Schultz quantum chain p2 |. These are chains 



with (m + n) L states invariant under the quantum superalgebra U q SU(m/n) p3| . 
In this paper, the (2/0), (1/1), (3/0) and (2/1) Perk-Schultz models will play a role 
with q real (\(3\ > 2) and the physical significance of the deformation parameter q will 
become apparent. We also give some new representations of the quotients. First we 
give non-hermitian representations (for q real), these are relevant for expressions like 
Eq. ( |1.9| ) and next, we give representations of the (2/0) and (1/1) quotients with 3 L 
states. Notice that in the last case the symmetries of the corresponding chain are not 
anymore U q SU(2/0) and U g SU (1/1). The knowledge of these symmetries is impor- 
tant because if two chains belong to the same quotient, their spectrum is in general 
the same but the degeneracies are fixed by the symmetries (a more detailed version 
of Appendix C is going to be published elsewhere). After finishing this appendix, the 
"mathematician" should go through Sections [||4] and have a close look at Section |5| 
where, in the simple case where the calculations can be done using free fermions, 
one illustrates the peculiar problematics of nonequilibrium statistical mechanics. It 
is stressed that (see Eq. (|1.S|)) the knowledge of the wave function (not only the spec- 
trum) plays a crucial role and since in the Bethe ansatz this knowledge is hard to get, 
the calculation of average quantities presents a new challenge. The "mathematician" 
should next have a look a Appendices A and B and skip Sections [7]-|9|. 

We sugggest to the "physicist" to read the paper in the chronological order. In 
Section |2| we consider quantum chains with L sites (we always take open chains). 
On each site we take a discrete variable f3 taking iV values ((3 = corresponds to a 
vacancy). We write the most general master equation describing bimolecular reac- 
tions. In Section [| we write down the corresponding one-dimensional Hamiltonian, 
see Eq. ( p..8| ) . A close related development connecting the master equation in discrete 
time to the transfer matrix formalism can be found in Refs. . 



Two-states models are considered in Section The elementary processes in the 
master equation describe besides diffusion, annihilation, coagulation and death pro- 
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cesses (these processes lead for large times to a state with vacancies only) also the 
reverse processes creation, decoagulation and birth. We first show that for pure 
diffusion processes which are left-right asymmetric, the Hamiltonian is just the q- 
deformed XXZ spin-1/2 Heisenberg chain. One discusses in detail the quantum 
chains corresponding to the different processes and one stresses the importance of 
the phase diagram of the XXZ Heisenberg chain in a Z magnetic field, especially the 
Pokrovsky-Talapov line. The learned reader will also notice the importance of non- 
hermitian representations of the braid group occuring in this type of problems. As 
is well known when a forward-backward process exists (annihilation and creation for 
example) through a similarity transformation, the Hamiltonian is hermitian. When 
all three forward-backward processes are allowed, this is not always possible. We 
derive the conditions on the rates in oder to get hermitian Hamiltonians. 

Section ^ illustrates in detail the simple example of annihilation only, with a rate 
equal to half the diffusion constant. This corresponds to the physical picture in which, 
when two particles are on neighboring sites, they always annihilate. As was already 
known, in this case all the calculations can be done using free fermions. Here we 
bring some new results. We consider the finite-size scaling of the problem (large time 
and lattice size with z = tL~ 2 fixed). We show that the finite-size scaling function 
exists (no logarithmic corrections). This result is important because if one accepts its 
general validity, it allows numerical estimates of critical exponents like in equilibrium 
problems. We also compute, for the first time, using the Hamiltonian formalism the 
density-density correlation function and stress the importance of the scaling limit 
r, t — > oo with u = r 2 /t fixed. 

In Section |5] we consider three-state models (two species of particles and vacancies) 
with Z 3 symmetry. Various integrable quantum chains occur which allows us to obtain 
some rigrous results. Nevertheless, as will be seen much work has still to be done. 
An interesting three-states model with Z 2 symmetry is mentioned in Appendix C. 
In Section |7] we derive and solve the condition on the rates in the master equation 
which gives chosen steady states. Based on these results, in Sections || we discuss the 
dynamics of the Ising and chiral Potts models. This takes us back to the quantum 
chains discussed in Sections |] and |6]. We look at the behaviour of the systems 
when the temperature is small (critical dynamics). This brings us to a problem of 
perturbation theory. Two cases occur. In the first, although the (spatial) correlation 
length diverges, the relaxation time stays finite. If we take the corresponding rate 
to zero and get critical dynamics, then the relaxation time depends on the (spatial) 
correlation length in a way which is independent of the remaining rates (universality). 

In Appendix A we study closer Eq. ( |1.9|) . We show how from the knowledge of 
the eigenvalues and eigenvectors of Hq one can compute the eigenvectors of H . In 
Appendix B we consider the example studied in Section |5| from a different point of 
view. We notice that the non-hermitian Hermiltonian corresponds to a representation 
of the Hecke algebra and using Baxterization, we derive the corresponding vertex 
model. It turns out that this is a seven- vertex model. This observation is relevant, 
since as shown in Appendix C (see also [25|) there are other examples of non-hermitian 
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chains (irreversible processes!) which satisfy the Hecke algebra and through this 
procedure one can find the wave functions using the Bethe ansatz. Section 10 closes 
the paper with some open questions. 



2 The master equation 

In order to write the master equation which describes a general lattice version of 
a reaction-diffusion process in one dimension, we take a chain with L sites and at 
each site i we take a variable $ taking N integer values = 0, 1, . . . , N — 1). By 
convention we attach the value $ = to a vacancy (inert state). We want to consider 
a master equation for the probability distribution P({j3}; t) with the following form 



dP{ffl;t) 



L-X 

<)f ^ 

Ul k=X 

N-X 



+ E' w e%(Pk, Pk+x)P(Px, ...,Pk + t,Pk+i + ™,..., /3 L] t) 



e,m=o 



(2.1) 



(k) 

where w er ^ are the transition rates and the prime in the second sum indicates that 
the pair £ = m = should be excluded. We assume all the additions performed on 
the Pi to be done modulo N. 

The advantage of this notation is that one can introduce discrete symmetries in a 
simple way. We shall assume hereafter that the system is homogenous which implies 
that the transitions rates are independent of k 

w { e k l(a, P) = w e , m {a, P) (2.2) 

for all k = 1, . . . , L — 1. The probability r^'^l that a state (7, 5) on two consecutive 
sites will change after an unit time into the state (a,P) is 

r%% = w,. a ^(a, P) ; (a, P) £ (7, 5) (2.3) 

The rates Wo,o(oc,P) are related to the probability that in the unit time the state 
(a,P) unchanges. From the conservation of probabilities, we have 

w 0fi (a,P) = Y^'w r , s ([a -r,P - s) (2.4) 

r,s 

where r = s = is again excluded. 

It is now trivial to check, using Eq. ( |2.3|) , that if for the iV-state model we want 
to have a Zn symmetry, only the functions 

w e>m (a,P) , £ + m = (modulo N) (2.5) 

will appear. This is the case of the annihilation model ( written &sA + A—»-® + $ 
where is a vacancy) which has Z 2 symmetry (N = 2) and for the two-species 
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annihilation reaction (N = 3 and the reaction is written as A + B — > + 0) which has 
Z 3 symmetry. In the former case, we assign to the vacancy the Z 2 quantum number 
and to the A molecule the number 1. In the latter case, the vacancy is denoted by 
0, the A by 1 and the B by 2. 

Parity conservation (left-right symmetry) is achieved if 

w e ,m(a, (3) = w m>e (f3, a) (2.6) 

Let us comment on the supplementary symmetries besides a possible parity invari- 
ance which exists in the four examples given in Section p]. The coagulation model has 
no symmetry at all. As we have seen, the annihilation model has a Z 2 symmetry. The 
two-species trapping reaction has a U(l) symmetry (since the number of A particles 
is conserved) and a Z 2 symmetry (give to the vacany and the state B a Z 2 -parity "+" 
and to the state A a ^-parity "-"). As mentioned before, the two-particle reaction 
has a Z3 symmetry. It has also a supplementary U(l) symmetry, since the difference 
of the numbers of A and B particles is conserved. 

It is often useful to make a change of variables in the master equation. In this 
paper, we define 

P(m;t) = mi3}M{i3};t) (2.7) 
where $({/5}) takes the special form 

mp}) = f[h {k) m (2.8) 

k=l 

From Eq. ( [2.1| ) ^ is a solution of the new master equation 

-w ,o(&,/3*+i)*({/3};*) 

+ £V$(ft, A+i Wi, • • • , Pk + *, Pk+i + m, . . . , p L ; t) 

t,m 

(2.9) 

with the pair £ = m = excluded and 

Wtlic 0) = yW( ;<^ m W °. 0) (2.10) 

where 

^ k \a,(3) = h^\a)h^ k+1 \f3) ; k = 1, 2, L - 1. (2.11) 

Notice that although the rates W£ >m (ct, (3) are link independent in general, the W^(a, 0) 
are not. The function $({/?}) has obviously to be non-zero and finite. Let us finally 
note that if we want all "molecules" to disappear at large times so that we are only 
left with vacancies, we must have the condition 

w e , m (-£,-m) = (2.12) 



UL fc=l 
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for £, m = 0, 1, . . . , N — 1. One can check that 



fe=i 



(2-13) 



satisfies indeed dP/dt = 0, with the convention that (3k = corresponds to a vacancy. 



3 Quantum chains corresponding to master equa- 
tions 



We shall now write the master equation (|2.9| ) in the form of a Schrodinger equation. 
In order to do so, on each site we define a basis in the space of N x N matrices E ki . 
The only non- vanishing matrix element of the matrix E kl is the one on the k th line 
and the I th column and this matrix element is equal to unity. Assuming homogeneity, 
the most general N L x N L Hamiltonian with only nearest-neighbor interactions can 
be written as 



L-l 



i=l 



where 



Hi =22 CL£ t . 



tm 



E r 



(3.1) 



(3.2) 



acts in the subspace <8> V^ +1 ^ of the A^ L -dimensional vector space. Besides the 
matrices E M , on each site it is also convenient to define the matrix F 



( 1 ••- \ 
1 ■■■0 



F 



(3.3) 



••• 1 

V i o o ••• o o j 

With these notations, the Hamiltonian density can be written as 

Hi = Ui — Ti 

where 

N-l N-l 

^ = Y! E W t %(a,(3) (E aa F e ) ® (e^F" 
e,m=o a,p=a 

N-l 

Ui = ]T wo,o(a,P)E aa ®E^ 

a,P=0 



(3.4) 

(3.5) 
(3.6) 
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Here is given in Eq. (|2.10|) and w 0i0 in Eq. (p.4|) . The Schrodinger equation 

replacing the master equation reads 

^\V) = -H\9) (3.7) 

The Schrodinger equation corresponding to the physical problem (without the simi- 
larity transformation Eq. ( |2.7|) ) is 

l\P) = -H\P) (3.8) 



H is obtained by taking w^ m (a } f3) instead of We :in (a, (3) in Eq. (|3.5| ). Obviously H 
and H have the same spectra but they are given by two different one- dimensional 
quantum chains which act in a N L Fock space (for N — 1 species) and their explicit 
form depends on the chemical reaction and diffusion process. We first discuss the 



Schrodinger equation (|3.8|) The ket state \P) is defined as follows p0| . Take an 
orthogonal basis in {/?} 

\{P}) = \Pi,...,Pl) ; <{/nil{/3}> = Vkw (3-9) 

then 

\P) = J2 p m;t)\{(3}) (3.10) 

The reaction-diffusion process is determined by the initial (t = 0) probability distri- 
bution P ({/?}) which defines the "initial" ket state 

\Po) = Y,p°m)m) (3.ii) 

{/3} 

The Hamiltonian H is in general non-hermitian and due to probability conservation 
(Eq. ( p.4|) ), it satisfies the remarkable relation: 

(0|F = (3.12) 

where the bra ground state (0| is 

<0|=£<{/3}l (3-13) 

{/3} 

From Eq. ( |3.12j ) it follows that the ground state energy is zero. Take now an observ- 
able X ({/?}) (for example the concentration of A particles in the coagulation model 
( |1.1|) ). Its average can be computed as follows 

<X>(t) = £X({/3})P({/3};t) 

= (0\X\P) = (0\Xe~ Ht \P ) (3.14) 
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Notice that that in nonequilibrium problems one studies the properties of the wave 
function which is already a probability and not quantum mechanical expectation 
values 

(0|X|0) (3.15) 

as one does in equilibrium problems. The thermodynamical (continuum) limit can 
be computed from Eq. fl3.14| ), taking the length L of the chain to infinity for a fixed 
time t. As discussed in detail in Section ||, a second limit (the finite-size scaling limit) 
is also interesting, where one takes both t and L large but keeps z = t/L 2 finite. 
If E\ and \^\) are the eigenvalues and eigenkets of H, we have from Eq. ( p.!4j ) 

<X>(t) = J2 a ^~ Ext (°l X I*a) (3.16) 

A 

where 

|Po> = E a A|^A> (3.17) 

A 

Thus the large t behaviour of < X > (t) is governed by the lowest excitations of H . 
If, instead of H we use H (see Eq. ( p.7|) ), then the averages have different expressions 



<X>(t) = J2X({P}M{P}M{Pht) = (0\X®e- m |tf > (3.18) 



where 



l*o) = E $ " 1 ({/ ? })^({/5})|{/3}) (3.19) 

There is another distinctive feature of nonequilibrium problems, as compared to 
equilibrium ones, and this is the concept of interaction range. As opposed to equilib- 
rium problems where the whole information is only contained in the Hamiltonian (it 
is nearest-neighbour interactions as can be seen from Eq. (|3.1|) ), in the nonequilib- 
rium case we have to give also Pq ({/?}). This probability can describe an uncorrelated 
homogenuous distribution like 

Po({P}) = f[fm (3-20) 

or a strongly correlated distribution when, for example, at t = the reactants are 
separated in space ||. The general properties of P({/3};i) and implicitly those of 
averages (like self-organization, critical dimensions or critical exponents) are going 
to be different. This can be understood when comparing Eqs. ( |3.14[ ) and ( |3.18| ). Let 



us assume that we give H and Pq describing a correlated distribution, we can make 
a similarity transformation to bring P to an uncorrelated distribution ty . After 
this transformation we will have to work with a new Hamiltonian H (in general with 
long-range interactions) with different physical properties. 



10 



4 Two-state Hamiltonians 

Since for two-states models instead of the basis E a, P one often prefers the basis of 
Pauli matrices, let us start by giving some useful identities 

E 01 g) E 10 + E 10 <g> E m = l -{a x ® a x + a v ® a v ) 

E 01 ®E 01 = ^(a x ®a x -a y ®a y + i(a x ®a y + a y ®a x )) 
E m ® E m + E 10 ® E 10 = ]^(a x ® a x - a y ® o y ) (4.1) 

We are now going to consider various Hamiltonians according to their symmetries 
and chemical properties. In the two-states models, we have A states and vacancies. 

4.1 Z 2 symmetric, parity non-invariant vacuum-driven pro- 
cesses 

By a vacuum-driven process, or inert-driven process, we mean reactions which end 
in a state with vacancies only. This implies from Eq. ( |2.12| ) 



w lt0 {l, 0) = u> O) i(0, 1) = Wl)1 (l, 1) = (4.2) 

This means that there is no production of pairs of A particles, in other words the 
process + 0^A + Ais forbidden. Since 

w lfi (a, (3) = w ,i(a, (3) = (4.3) 

because of the Z2 symmetry, we are left only with the processes 

a. ) annihilation, with the rate wi^O, 0) 

A + A -> + (4.4) 

b. ) diffusion to the right, with the rate u>i,i(0, 1) 

A + dS^dl + A (4.5) 

c. ) diffusion to the left, with the rate u>i,i(l, 0) 

+ A->A + (4.6) 



We use now Eqs. (|3.4|, 3.5, |3.6|), take into account that 



E 00 F = E 01 , E n F = E 10 , a z = E 00 - E n , (4.7) 
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and find 



T 



¥><*>(!, 1) 



w 1)1 (0,0)E 01 ®E 01 



<^)(0,0) 

, y (i) (Q,i) 
V } (i,o) 



Ui 



- (wi,x(0, 1) + w 1A (l, 0) + iu^O, 0)) 1 ® 1 
+^ (^(0, 0) - w ltl (0, 1) - 0)) <t* ® (7 2 



101,1(0, 0) (a 2 ® 1 + 1 ® a z ) 



+ - K,i(l, 0) - 101,1(0, 1)) (a 2 ® 1 - 1 ® a z ) 



(4.8) 



Up to this point we still have freedom in the choice of the function ^ k \(3k-> Pk+i) in 
Eq. Q2.8Q . We now take 



¥>«(!, 0) 



<^)(0,1) \ u;i,i(0,l) 



^i,i(l,0) 



q ; fc = 1,2, ...,L - 1, 



(4.9) 



and it is clear that g is real since the rates are real. From Eqs. ( |4.9j ) and ( |2.11| ) we 

get 



^ fc )(l) 



/*( fc )(o) 

where A is an arbitrary parameter. The diffusion constant D, expressed as 

D = y/w ltl (0, 1)^(1,0), 
fixes the time scale. We then finally obtain 

H = D(H + Hi) 

where 

#o = [of of+i + ofof+i + Atrf <r? +1 + (1 - A') (of + of +1 ) 



(4.10) 



(4.11) 



(4.12) 



i=l 



q-q 



<Ja — a. 



i+l 



2 A' - A - 2 



L-l 



^1,1 (0,0) x I —2/ / -01 ,'01 



\ 2 D 

q + <i 



i+l 



i=l 
-1 



wi.i(0,0) A / = 1 _ ^i,i(0,0) 



2D 



2D 



(4.13) 
(4.14) 

(4.15) 
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Let us now discuss the structure of this Hamiltonian. H is hermitian, £7(1) symmetric 
and its properties are going to be discussed shortly. Hi is non-hermitian and has only 
a Z 2 symmetry, corresponding to the transformations 

a? -> -af , aV^-aV , 0*^0* (4.16) 

which makes H only Z 2 -symmetric (as expected) and non-hermitian. 

Let us discuss the case of pure diffusion, that is Wi 5 i(0,0) = 0. We have Hi = 
and 



■ ' ' ' ' .,</<' 1 9 + 9 1 „z„z '/ '/ ' (■, _z \ 'I ■ 'I 



'x 1 ^.y ^.y j. • j. „z 1 ~z ^ 

o i a i+1 + of of +1 + — - — a i a i+1 — [a t - a. 



't+i 



1 L-l 

(4.17) 

This Hamiltonian is the U q SU (2) symmetric Hamiltonian of Pasquier and Saleur p6| 
where the deformation parameter q has a clear physical interpretation, see Eq. ( |4.9| ) . 
Moreover, as shown in Appendix C, Hq is the sum of generators of a given quotient of 
a Hecke algebra (called a Temperley-Lieb algebra). As discussed in detail in Appendix 
C, for any number of species, asymmetric diffusion defines various quotients of Hecke 
algebras (we think this is an important observation). As a corollary one can show 
that for any number of species asymmetric diffusion is always massive. 

We now consider the case 1^1(0, 0) 7^ 0, that is the full nonhermitian Hamiltonian 
Eq. ( |4. 13| ) which has a special property. From the expression ( f4.15| ) of Hi we see 



that one can take the factor A 2 arbitrary without changing the spectrum. This is 
certainly so, since this factor only serves to parametrize a similarity transformation. 
On the other hand, one can also see by direct inspection of the matrix elements of H 
that the characteristic polynomial of H is not affected by the existence of H\. Let us 
briefly present the argument since we are going to use this repeatedly later on. If C 
and D are square matrices and if A is a rectangular matrix, it is well known that 

detCdetD (4.18) 

Now, Hq has a (7(1) symmetry and has thus a block- diagonal form. The corre- 
sponding quantum number labelling the blocks is the number of A particles present. 
Acting with Hi reduces this quantum number by two and this plays the role of X in 




Eq. ( |4.18| ). Thus the above determinant formula can be applied and the independence 
of the spectrum of H from the operator Hi follows. 

Let us choose ^1^(0,0) such that A = 0. This corresponds accidentally to the 
choice of all the previous numerical simulation studies of the one-species annihilation 
process where one took q = 1 and Wi i(0,0) = 2D. In this case H can be studied 
in terms of free fermions. This will be done in Section |5| where some new results 
are presented. In particular we find that for q — 1, Hq can be written in terms of 
fermionic number operators as follows 

H o~Y,(t) a Uk (4.19) 
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which describes a Pokrovsky-Talapov phase transition p7| . Notice the quadratic 
dispersion relation which is expected in any diffusion problem. 

In Appendix B we consider a complementary point of view which is related to the 
problem of integrability of quantum chains with the structure Eq. ( (4.12j ) in which Hq 
is integrable (in the present case it corresponds to a six-vertex model) and Hi does 
not affect the spectrum of H. We show that H is related to a seven- vertex model for 
which the Yang-Baxter equations are valid and hence H is integrable. 

We now consider the case A / 0. This brings us, see Eq. ft4.13|) , to discuss 
properties of the Hamiltonian 

H' = -\ E hX+i + *Mn + A «i + h {°l + <fi)l ( 4 - 2 °) 



This Hamiltonian is integrable G3J and its phase diagram is shown in Fig. 1. For 



h = 0, the system is massive with a ferromagnetic ground state if A > 1, massless and 
conformally invariant if — 1 < A < 1 and again massive with an antiferromagnetic 
ground state if A < —1. For a given h, the system is massive ferromagnetic if 
A > 1 — h, then undergoes a Pokrovsky-Talapov transition at Ap^ = 1 — h, is in a 
massless incommensurate phase for A c < A < 1 — h and reaches again the massive 
antiferromagnetic phase if A < A c . 

It can be shown |29| that in the continuum, along the Pokrovsky-Talapov line, 



the spectrum of H' is, up to normalisation, given by Eq. ( |4.19| ) for any A. In other 
words, the system is massless with a quadratic dispersion in the momentum k. 

With this knowledge in hand, let us discuss some properties of H . Since A' < A, 
we have h > 1 — A and the system is massive. If however, q — 1, we have A' = A 
and we are on the Pokrovsky-Talapov line where the system is massless. There are 
some immediate questions to ask about this system. If we define 

«-a^-i (4.2i) 

and if r denotes the relaxation time, we are interested in the exponent rj 

t ~ e~" (4.22) 

It can be shown that rj = 1 is independent of A using the standard lore of Hecke 
algebras, see Appendix C. This result is to be expected since we can change the value 
of A by changing D, see Eq. ( |4.15|) . However, changing D merely changes the time 



scale. More generally, the full finite-size scaling form should read 

r = L m F(x) , x = eL m (4.23) 

in the simultaneous limit e — > 0, L — > oo with x kept fixed. Concerning the concen- 
tration of A particles 

ca = j E El 1 (4.24) 
L i=i 
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in the critical regime (q = 1), we are interested in 



<c A > (t) = L x c 



tlr 



(4.25) 



for both t and L being large with z fixed. The exponent 2 in the definition of z 
was taken because of the quadratic dispersion relation Eq. ( f4.19|) . In Section || we 
shall find x = — 1. Similarly, correlation functions can be introduced and calculated. 
We shall present the analysis in Section for the case A = 0. For the case A ^ 0, 
although we expect on physical grounds the same results, the explicit proof is still 
missing. The Bethe ansatz equations for the wave equations are known only for H 
but not for H. In Appendix A it is shown how in principle the knowledge of the 
eigenfunctions of H can help to find the eigenfunctions of H . In a different approach 
one could start with the seven- vertex model (which is not limited to the case A = 0), 
and perform the Bethe ansatz there. The wave functions thus found could be used 
for H. Whether this whole program is manageable remains to be seen. 



4.2 Z2 non-invariant, parity invariant vacuum- driven pro- 
cesses 

From now on, we take always the left-right symmetric case, see Eq. ( p.6|) , in all 
reaction rates, in particular wi,i(0, 1) = u>i,i(l, 0) = D. We choose units of time such 
that D — 1. To the processes studied before we now add the following 



d. ) coagulation, with rate u>i,o(0, 1) 

A + A -> + A 

e. ) decoagulation, with rate u>i,o(l, 1) 



+A^A+A 



f.) death, with rate iOio(0, 0) 



We take <p®(0, 1) = (pM(l, 0)(z = 1, 2, L—l) and get 

^(0(1,1 



(4.26) 

(4.27) 
(4.28) 



T; 



E 01 ®E W + E W ®E 01 ) + 



w hl (0,0)E 01 ® E m 





(1,0) 




(0,0) 


K 


(1,1) 




(1,0) 


K 


(1,0) 




(1,1) 



¥>W(0,0) 

' :Wl , (0,0) {E 01 ®E 00 + E 00 ®E 01 
^> li0 (0, 1) (E 01 ® E 11 + E u ® E 01 



^i,o(l,l) (E w ® E 11 + E u ® E w ) 



(4.29) 
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Ui 



1 ft n\ , 1 fi i\ i ( ™o,o(M) 
-w ,o(l, 0) + -w 0> o(l, 1) + I 

u,uv ; 1 ® (j 2 + a 2 ® 1 

4 



^o,o(l,0) > 



CJ Z <g> <J Z 



(4.30) 



where 



w ,o(l,l) = 2™ 1) o(0,l) + w 1 , 1 (0,0) 
w ,o(l,0) = wi, (0,0)+w 1 ,o(l,l) + l 



(4.31) 



Thus the Hamiltonian depends on four parameters. Two cases have to be considered 
separately. In the first case we have only coagulation, in the second case we have 
coagulation and decoagulation. If we have no decoagulation, n»io(l, 1) = 0. Then 



H = H + H 1 



(4.32) 



where 
#o = 

Hi = 



\ E [«H + «Mu + ^M+i + (1 - A') 0* + + 2A' - A - 2 

(4.33) 



i=l 



L-l 

E 

i=l 



¥»W(0,0) 



™ 1 , 1 (o,o)ifi^ 1 + 



y (0 (i,Q) , 

y»W(0,0)' 



-.00 7-.01 



+ ^^^(0, 1) + E^l, 



(4.34) 



A = 1 - «7i >0 (0, 1) - |wi t i(0, 0) (4.35) 
A = A' + U7 li0 (0,0) (4.36) 

We have checked using the same arguments as above that the spectrum of H is 
independent of the presence of Hi. Since A' < A, we are in a massive phase. The 
system is massless, i.e. A = A', if wi )0 (0, 0) = 0, corresponding to the absence of 
death processes. A straightforward calculation shows that the relaxation time r scales 
like 

r-M" (4.37) 
This can be seen by repeating the same arguments as for Eq. ([4.22 ). The result 



Eq. ( 4.37 ) can be understood easily if one keeps in mind the time evolution of the 



system. If one considers the later stages of the process when few "molecules" are left 
few of them will meet and annihilate. On the other hand, any one of them can "die' 
since this is an individual process. 
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In the absence of the death process some results are already known for the pure 
annihilation process [12[ or the coagulation process ITT] and it remains to be seen 



how much more can be learnt using the quantum chain formulation. 

We now consider the case when the coagulation and the decoagulation processes 
coexist and we have all terms in Eq. (|4.29| ). Then little can be said. One interesting 



case 



30| is when apart from just diffusion, only coagulation and decoagulation are 



present, that is w^i(0,0) = u>i )0 (0,0) = 0. We choose 



¥>(*>(!, 0) \ ^,0(0,1 



u>i,o(l, !) 



1,2,...,L-1). 



(4.38) 



We stress that this transformation is singular when the decoagulation rate wx t o(l, 1) 
goes to zero. We get a hermitian Hamiltonian H = H + H\ with H being given by 
Eq. dg|) with 



A' = 1-^0(0,1) 
A = A' + w lt0 (l,l) 



and 



L-l 



H 1 = -Y J >/(A-A')(l-A') {at Ell, + V+i) 



(4.39) 
(4.40) 



(4.41) 



where A' < A and again we suspect a massive phase, at least when A' is close to A. 
Note that now Hi is hermitian and its presence does change the spectrum of H. This 
change is in fact quite important. If we look directly at the master equation, besides 
the trivial stationary probability distribution Eq. ( |2.13|) there is a second one which 
satisfies dtP = and is given by 



P({PY,t) = Up(Pk) 



fc=i 



where 



cxp 



(M- 



In 



wi,o(o,i; 
^1,0(1,1; 



-[L - In 1 + 



wi,o(l, 1) 
wi,o(0,l). 



(4.42) 



(4.43) 



The probability distribution Eq. ( |4.42j ) corresponds to a one- dimensional Ising model 
defined on the dual lattice and the site variable (3k corresponds to a link variable of 
the dual lattice, see also Section ||. 

The existence of two stationary probability distributions satisfying d t P = makes 
the Hamiltonian H have a degenerate vacuum. The existence of a decoagulation 
process, no matter how small it is, implies, if the annihilation process is absent, 
that the stationary configuration should correspond to Eq. ( 4.42| ), unless the initial 
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state is the empty lattice. If we now turn on the annihilation mechanism, the final 
configuration will be the empty lattice. Since we always have massive behaviour this 
transition between stationary states will correspond to a first-order nonequilibrium 
phase transition which may be accounted by the fact that no symmetry of H is 
broken. This can now be studied by performing perturbative calculations in the limit 
of a small decoagulation rate. 

We are aware of the following fascinating puzzle. Although the Hamiltonian 
describing the coagulation and decoagulation processes is not known to be integrable 
it was shown in that the gap-probability function for this model can be computed 
exactly and this is an indication that the model is integrable. This is at least the case if 
one considers the creation-annihilation model with diffusion which will be discussed 
next (recall that Glauber's solution was discovered without using the Jordan- 
Wigner transformation). We will return to this problem in a future publication (see 
also Section |]). 



4.3 Z2 and parity invariant, not vacuum-driven processes 

The full master equation without any restrictions includes besides the processes con- 
sidered so far the following two, see Eq. ( |4.2| ) 



g.) creation, with rate i0i,i(l, 1) 



+ 



A + A 



(4.44) 



h.) birth, with rate i0i,o(l, 0) 



From probability conservation, we have 



A + 



wo, (0,0) = u>i,i(1,1) + 2wi j0 (1,0) 
w ,o(l,0) = u;i,o(0,0)+twi i o(l,l) + wi,i(l,0) 
u; ,o(l,l) = 2u;i i0 (0, 1) + u;i,i(0, 0) 



The Hamitonian is 
H 



L-l 

E 

i=i 

1 



- Ko(0, 0) + 2«j ,o(1, 0) + wo,o(l, 1)) 



+ - Ko(0, 0) + i0 O ,o(l, 1) - 2i0 O)O (l, 0)) <t?0? +1 



+i(«, 0i0 (o,o)- 100,0(1, i))(^ + ^ +1 ) 

-1^(1,0) (E^ + EfE^) 



(4.45) 



(4.46) 
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-A-V,o(l,0) (e?E% 1 + EFe£ 1 
-A«;i,o(0,1) (E^Ell. + El'E: 
-A-V,o(l,l) (EFE^+EfE, 



10 

i+l 



(4.47) 



where we have choosen 



A, (k = 1,2, ...,L) in Eq. ( |2.11| ). This Hamiltonian 
has an interesting property. It can always be brought to a hermitian form through 
a similarity transformation if one out of the three possible forward-backward reac- 
tions (annihilation-creation, death-birth, coagulation-decoagulation) take place. If, 
however, the following relation between the rates is satisfied 



101,1(1,1) /Vl,o(M) 



wi,i(0,0) Ui.oM 



Vq(1,0) ' 
wi, (0,0) ( 



(4.48) 



H can be made hermitian even when all three pairs of reactions occur. The physical 
significance of Eq. (|4.48|) has still to be explored. 

We shall confine ourselves to the Z 2 preserving processes which include besides 
diffusion only annihilation and creation processes. We make the transformation 



«(1,1) 



^,i(l,l) 



(k — 1, 2, ...,L — 1), 



^( fc )(0,0) ^i,i(0,0) 
and find H/w\ t i(l, 0) = H + H% where 

#o = -5Ek^i+^M + i + Aa^ +1 + (l-A')(^ + a, 
Hi = ~\ £ V / ( A ' - A)(2 - A - A') (ofof +1 - o\oU 



(4.49) 



t+i 



A - 2 



(4.50) 



A' 
A 



«»i i(0,0) -^i,i(l,l) 



2wi,i(l,0) 



A'- 



w M (l,0) 



(4.51) 
(4.52) 



This is just the Hamiltonian corresponding to the kinetic Ising model with purely re- 
laxational dynamics |3T] and H is hermitian. The limit when the elements of H% are 
small does need some care. If the annihilation and creation rates are equal, A' = 1 
which allows us to do perturbation theory in the coupling 1 — A of Eq. ( |4.50| ). How- 
ever, if the rates are not equal and one just takes iui i(l, 1) to zero, the limit for the 
eigenvectors is singular as is the transformation Eq. ( 4.49Q , although the eigenvalue 
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spectrum is not affected. Therefore, the heuristically appealing relationship [32] be- 
tween the kinetic Ising model with Glauber dynamics and the annihilation process 
should be considered with care if one works (like in Ref. |3lfl ) with the hermitian 
formulation given by Eq. ( }4.50| ). 



5 Exact solution of a two-state system 

In this section, we shall discuss the dynamics of the two-states system described 
by the Hamiltonian Eq. ( |4.12| ) in the free fermion case A = 0. (For a connexion 
between this Hamiltonian and a seven-vertex model, see Appendix B.) This model 
describes A particles diffusing and annihilating in pairs. The motivation for treating 
this particular case in great detail is twofold. 

Firstly, we would like to show that a finite-size scaling theory can be formulated for 
nonequilibrium systems. This point is particularly important if (in contradistinction 
to this example) the problem is not exactly solvable. In that case, one will have to 
solve the problem numerically for finite chains and then extrapolate the results to 
infinite systems, as one does for equilibrium problems [|33| (see also Appendix A). 
One expects that the finite-size scaling proven for this model is of general validity. 

For instance, the concentration of particles per site <n > / L obeys the relation 

c=L x <j>{z) (5.1) 

where z = 4Dt/L 2 and 4>(z) ~ z~ a as z — > 0. Moreover, x is related to the large time 
behaviour of c ~ t~ a , namely x = —2a. 

Secondly, we want to stress the interest of the (connected) two-point function 

G(r,t) =< n m n m+r > - < n m >< n m+r > (5.2) 

where m indicates the position of a lattice site. Here, we shall pose two different 
questions about G(r,t). First, we ask for the large time behaviour (t — > oo) if r is 
kept fixed and look for the exponent v r 

G{r,t)~f r (5.3) 

Secondly, we consider the scaling limit where simultaneously r, t — > oo such that 
u = r 2 /t is fixed and define a critical exponent y 

G(r,t) ~ r~ y g(u) (5.4) 

where g is a scaling function. We believe that for nonequilibrium problems the 
exponents ?/'s might turn out to be as fundamental as the critical exponents of two- 
dimensional equilibrium statistical mechanics, since the information given by the 
initial probability function should be hidden in the scaling function g(u) and conse- 



quently the y's should be universal (see also Ref. |32|] ) 
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Finally, we would like to present powerful techniques which allow us to reduce 
the calculation of nonequilibrium averages, see eq. ( |3.18| ), to vacuum expectation 



values, thereby extending the approach initiated by Lushnikov to the calculation 
of correlation functions. 

The quantum Hamiltonian H ( |4.12j ) contains the parameter 



\ (q + q' 1 ) - 1 > (5.5) 



£ 2 

where e measures the left-right asymmetry of the diffusion process. In the remainder 
of this section we take e = 0. We do not make use of the change of variables given 
by the function (3) described in the last section. In terms of Pauli matrices 

a ± = (a x ± ia y ) /2, the Hamiltonian is equivalent to 

L 

H =~ D Y. K>n+1 + °n°n+l ~ 2 °n°n + 2 °n°n+l] ( 5 - 6 ) 
n=l 

where D is the diffusion constant. In this section for simplicity we only use periodic 
boundary conditions. 

5.1 Diagonalization and the generating function 

The diagonalization H is fairly standard. By restricting ourselves to the case L = 2M 
even we define fermionic variables via the Jordan- Wigner transformation 

a+ = a + (-l)^™<» a " flm (5.7) 
as well as their Fourier transforms 



L -l/2 e -in/4J2 e - 

q 

L -l/2 e +w/4£ e -i,n a + ( 5 .{ 



1 

1 



In terms of these fermionic variables the Hamiltonian Eq. ( |5.6| ) takes the form 
H = —2D ^2 [(cos q — 1) (a^a q + at q a^^j + 2 sing a q a^ q 

0<q<ir 

- 5 g)7r 2a+a n .] (5.9) 
where q takes the values 

2k — 1 

geven = 2M n , k = 1, 2, . . . , M 
2k 

g ° dd = 2M n ' fc = ' 1 '--- M ( 5 - 10 ) 
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if the total number of fermions Af = J2m=i a m a m is an even or odd number, respec- 
tively. The Hilbert space associated to Eq. ( |5.9| ) can therefore be splitted into two 
disjoints sectors depending on the parity of the fermion number. The sector with M 
odd (TV even) will be related to the case where the ground state has a single particle 
(no particle). 

Looking again at Eq. ( [B~9f) , we note that H = J2o< q <wH q has block- diagonal 
structure. The time-evolution equation for wave functions d t $> = —HQ can be broken 
up into separate equations for each block 



m = n , ^ 

0<5<7T 



-H q $ g ; (0 < g < tt). 



(5.11) 



Each block is generated by acting with a q and a± q on the state \vac) without fermions, 
i.e. a q \vac) = 0. If q ^ 0, tt, H q can be written as a 4 x 4 matrix. In the basis 



|0) = \vac) , |1) 



| vac) , 1 2) = at | vac) and |3) = a+al \vac) it takes the form 



-2D 



( 



V 



cosg 



cosg — 1 



—2 sin q \ 



2 cos q — 2 J 



(5.12) 



while for q = 0, 7T, the blocks H and are already diagonal. The general solution 
of Eq. (|5.11| ) corresponds to 



%{t) = (a q {t)a+at q + -y 9 (t)a+ + 5 q {t)a± q + (3 q {t)) \vac) 
where the functions a, f3, 7, 5 satisfy the differential equations 



(5.13) 



a q (t) 



AD (cosg — 1) a q (t) 
—AD smqa q {t) 
2D (cos q -l)j q (t) 
2D (cos q - l)5 q (t) 



(5.14) 



The solutions of these equations are promptly derived 



a q (t) 



a q (0) exp [-ADt (1 - cos q)} 

/3,(0) + cot| {a q (t) - a q (0)} 

7 9 (0)exp [-2Dt (1 -cosg)] 
S q {0)exp [-2Dt{l - cosg)] 



(5.15) 



For an arbitrary initial condition the general solution will be given by linear combina- 
tions of the wave functions Eq. ( |5.11| ). A general study considering arbitrary general 
initial conditions is not straightforward and is not in our intentions. Since our main 
interest here is the finite-size scaling behaviour of the system, we will hereafter choose 
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the simple initial condition where we have no vacancies present. This condition was 
used by Lushnikov |34j] and corresponds to the situations where we have a single wave 
function (gUp , i.e. ip(t) = $(t) with 

a,(0) = -1 , (3 q (0) = 7,(0) = 5 q (0) = 0, (5.16) 

where q takes the values q even m the set ( |5.10| ), since in this case the fermion number 
M is even. 

In order to calculate the concentration of particles and correlations, following the 
work of Lushnikov [3^] it is convenient to introduce generating functions. Since (see 
Section H) the wave function takes the role of the probability distribution introduced 
earlier, a generating function convenient for our purpose is defined by 

F({z}, t) = F( Zl , ...,z L ,t)= E{c} ^ {c} . t) (5-17) 

where n, = 0, 1 and Zi denote the occupation number and fugacity at the site i, C 
is a configuration {rii, . . . , n L } of occupied or empty sites. In Eq. (|5.17|) we take 
into account that ip does not, in general, satisfy the normalization condition of a 
probability distribution. Using the results of Section || in particular (3.14), we can 
bring the calculation of F({z},t) to an equilibrium problem, i.e. 

F({z},t) = Fo(*)En«({C};t) (5.18) 

{C} i=l 

= F (t) (vac\ex P (j2z £ a~(£)) ]J (a q (t)~a + q ~a + _ q + (3 q (t)) \vac) (5.19) 

where we have already used the initial condition Eq. Q5.161 ) and F (t) is determined 
from the normalization condition F(l, . . . , 1, t) — 1 for all t. In order to understand 
Eq. ( |5.19| ), it is useful to consider the identity 

exp I z m a~(m) J = Yl z mx o~ im x ) ■ ■ ■ z mk a~~(m k ) (5.20) 

\m=l / k=0m\>m2>...>m k 

which we can derive by expanding the exponentials and using the relations (cr~(m)) 2 = 
and [er~(m), cr~(m')] = 0. It is now clear that (vac\ exp fe£=i zi&~(£j) corresponds 
to the state (0| of Section || 

In order to calculate F({z},t) we observe that due to the ordering m x > m 2 > 
. . . > m k we get 

(vac\ a~ (mi)o~ (m 2 ) ■ ..a~(m k ) = (vac\ a mi a m2 ...a mh (5.21) 

On the other hand, we also have from Eq. 

2 



at a 



q -q t 

m n>m 



E E sin(g(n - m)) a+a+ (5.22) 
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Because of the product structure Eqs. (|5.11| , |5.13| ) of 0, the generating function can 
also be written as F = II g >o-^<r Each factor contains only two fermionic creation 
operators at most. So we only need to calculate 

2J E ( vac \ cr~ (rn)a~ (m')d^ a+ \vac) 

ml m>m' 

= T E E E [sin(g(n-n')) (vac\a m a m ,a+a+, \vac) 

n' ,m' m>m' n>n' 

2 

= y E E E [ sin (<?(™ - n ')) ($n,m>5n>,m ~ S n ,m8n',m')} 
n' ,m' m>m' n>n' 

= -|E E 8in(g(n-n')) (5.23) 

n' n>n' 

In the sequel, we shall need two identities, which are obtained using Eq. ( |5.10|) , with 
m fixed 



sin(g(n — m)) 

i>m 

sin(g(m — n)) 



q cos(g(m — 1/2)) 

COt 1 ; — 

2 sin(g/2) 
q cos(g(m — 1/2)) 



cot 



2 sin(g/2) 
We now combine Eqs. fl5.18| ) - ( |5.24| ) and get 



F{{z},t) = F (t) II KW •f-yEE 2n«mBin(g(7i-m))J + 



(5.24) 



(5.25) 



We fix .Fo(t) from the condition -F(l, . . . , l,i) = 1 for all t. The final result for the 
generating function is 



_ tt ( ' (zi ^ m ^ n>m ZnZm sin ( q ( n ~ m ))) + 

(l4,t) "iU ^ -a,(t) cot § + 



(5.26) 



5.2 Calculation of the mean concentration 

The mean number of particles at an arbitrary site m for a chain of length L is now 
obtained from 



d 

< n m >= z m - — F({z},t) 

OZm, 



, ^ / -2a q (t) cot § 
, ZL=1 ~ 9 >o V ^t § + 



(5.27) 



Due to the translational invariance of the particular initial state considered here 
(Eq. ( |5.16| )) the concentration of particles per site is simply 

c =< n m > (5.28) 

Substituting Eqs. ( |5.13| , [5.16Q into (|5.27|) we obtain 

c(t) = 2L- 1 exp (-4Dt (1 - cos q)) (5.29) 

<j>0 
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5.3 Large-time behaviour and finite-size scaling 

We now discuss the late stages dependence of the concentration c(t). In the limit 
L — > oo for t fixed 



c (t) = - f e- iDt ^- cosg) dq = e- 4Dt I (ADt) (5.30) 

7T JO 

where Iq(x) is a modified Bessel function. For Dt ^> 1 one gets 

c(t) ~ (8nDt)- 1/2 (5.31) 

Let us now analyse the regime of finite- size scaling where L — > oo, t — ► oo but 
tL -2 fixed. In this case we obtain 

2 f:^ V (-2^z(k- l -]\ (5.32) 



-e 2 (0,27ixz) (5.33) 



= |^(l + 2f(-l)'exp(-£V(2,))) (5.34) 

where 02(z, t) is a Jacobi theta function and 

z = ADtL- 2 (5.35) 

is the finite-size scaling variable. Eq. ( |5.34 ) is obtained from ( p. 32] ) by using the 
Poisson resummation formula. From Eqs. ( |5.32| - |5~3"4| ) and fl5.1| ) we get the critical 
exponent x — — 1. 



5.4 Fluctuations around the mean values 

Next, we briefly discuss the fluctuations around < n m >. For simplicity, we take 
z\ = . . . = zl = z in Eq. ( |5.17 ). Then the generating function F = F(z,t) and we 
have 



<n m > 



2 = 1 



<n m > 2 +2 <n m > /L- —J2 



-2a q {t) cotf 



L &o V cot ! + 



(5.36) 



We get 



< n m > - < n m > 



2/L (e- 4D '/ (4Dt) - e- 8Ut I {8Dt) 



-8Dt , 



1 1 



Ly/2M V v^y 



(5.37) 
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as t — > oo. This means that for fixed, but large t 

<5n 2 m > / <n m > 2 ~ t l / 2 L- 1 (5.38) 

which accounts for the fact that as the time grows the annihilation of couples of 
particles induces larger fluctuations in the particle concentration. 



5.5 Correlation functions 

We now turn towards the calculation of correlation functions. As an example which 
illustrates the general technique, we calculate the connected correlation function de- 
fined in Eq. fl5.2|) . From the definition of the generating function (|5.17 ) we can write 



d 2 

G(r,t) =< n m n m+r > - < n m > 2 = - — In F({z},t) 

OZ m OZ m -\- r 

A straightforward calculation gives us 



(5.39) 

z 1 = ...=z L =l 



G(rt) = W( -^W*) )-±y( - a ^) cot l ) (5 40) 
[,) L^ [-a q (t) cot * + (3 q (t)) L^ Q {-a q (t) cot I + (3 q (t) J K ' 1 



which does not depend on m due to the fact that the Hamiltonian Eq. (5.6) as well 
as the initial probability distribution are translational invariant. 

Analysing the relative importance of the two terms contributing to G(r, t), we see 
that in the limit L — > oo, the first term is of order unity, while the second one is of 
order 0(L~ l ). 

This correlation can be evaluated in the same way as we did for the concentration. 
Let us consider initially the case of the correlation of next-neighbour particles. We 
find 

G(l,t) = j - cos g)exp[-4Dt(l -cos q)) 

L q>0 

= e- iDt [I {ADt)-h{ADt)] 

~ n{%TtDt)-* /2 (5.41) 

as t — > oo. The same result can be obtained using Glauber's dynamics for the Ising 
model at T = (see also Section and considering domain walls two-point functions 
|3"5| . The initial condition we are considering in the diffusion-annihilation problem 
(no vacancies) corresponds to a fully ordered antiferromagnetic state in the kinetic 
Ising model. On the other hand, for arbitrary but finite r, we obtain analogously 

G(r,t) ~rG(l,t) (5.42) 

This result reflects the fact that at late times it is more difficult to find particles close 
to each other than far apart. From Eq. (|5.3|) this gives us the exponents v r = —3/2. 
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Finally, in connection with Eq. (|5.4| ), let us consider the limit where r — > oo, 
t —>■ oo but u = r 2 /t stays finite. This leads to 

G ~ 7r(87r J D)- 3 / 2 r- 2 M 3 / 2 e- u /( 2D ) (5.43) 

and we read off the exponent y = —2 and recognize the scaling function g(u) = 
71(811 D)~ 3 ^ 2 u 3 ^ 2 e~ u ^ 2D \ see Eq. ( |5.4j ). We remark that if we consider the domain 
walls two-point function in the zero-temperature kinetic Ising model with Glauber 
dynamics for an initial probability distribution with zero magnetization (which 



is a different initial condition than that we are considering here), we obtain the same 
result. This is an indication in favour of our conjecture that the y exponents are 
universal. 

6 Three-state models with Z% symmetry and va- 
cuum-driven processes 

As noticed in Section |], behind many processes in the two-states models there is the 
XXZ model in a magnetic field. We shall now show that there are several integrable 
quantum chains which play a similar role for three-states models with Z% symmetry. 
We shall restrict ourselves to the Z 3 -symmetric case for simplicity. This already 
contains the two-species annihilation process but not the Z 2 symmetric trapping 
reaction to which we hope to return in a future publication. Let stand for a 
vacancy, 1 for the particle A and 2 for the particle B. 

We use the general results of Section || and the conditions Eqs. ( ^75l , |2.12|) and 
derive the Hamiltonian 

H = U-T (6.1) 



T 



U 



L-1 
i=l 

+ W2 , 1 (i,o)A; Jil C: 1 

+w 2 ,i (o,: 
+w 2 ,i(i, 

+W 2 ,1 (0, 

+w 2 ,i{2,: 

+ w 2 , 1 (l,2)E i ul E% 1 

L-1 

[Ki(2, 0) + w 1>2 (0, 2)) EfE] 1 ^ + (w 2>1 (0, 1) 

i=i 



2)E? 2 E™ 1 
I )^°^i 

21 zpOl 

ii^" ,2 4ii 

2)Ef 1 E^ 1 

llO 77.20 



y i+l 

+ w h2 (0, 
+ w 1>2 (l, 
+ w 1>2 (2, 
+ w h2 {0, 

+ Wi >2 (l, 

+ w lt2 (2, 
+ ^1,2(1, 
+ w lt2 (2 



1) ^1 

m^E 2 ^ 

2) £fi^i 

^ l)^ 2 °^i 



(6.2) 



^(1,0))^^ 
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+ (w 2A (l, 0) + w lt2 (2, 2)) E™EH x + (t« 3)1 (2, 2) + u, 1>2 (0, 1)) E^E™, 
+ (tw 2 ,i(l, 1) + Wl , 2 (2, 0)) E°°E% X + (tw 2)1 (0, 2) + w 1)2 {l, 1)) Ef 2 E™! 
+ Ki(2, 1) + w 1)2 (0, 0)) Ej^E^ + (tw 2)1 (0, 0) + iu lf2 (l, 2)) EfE 1 ^ 
+ (w 2 , 1 (l,2) + w 1>2 (2,l))E™E™ 1 . 

Notice that for the time being we do not use any equivalence transformation 
we did in Sec.(p. 

We first consider the case of asymmetric diffusion 



A + B 
5 + 
A + (H 
B + A 
+ 5 



and let us assume that we have 



A + with rate tw 2 ,i(l,0) 

B + A with rate w 2 ,i(2, 1) 

+ 5 with rate w 2 ,i (0,2) 

+ A with rate iui >2 (0, 1) 

A + B with rate w 1)2 (l,2) 

5 + with rate w x , 2 (2,0) 



(6.3) 
TSh as 



(6.4) 



By defining 

we now obtain 
H = 



™ 2l i(l,0) 
^'1,2(0, 1) 




-1 



q + q 



' i± f^ E KE 

Z a=0 



w 1}2 {2,0)=w 2)1 {2,l) 




9 = \ TT- 



b>a 



b<a 



aa 
+ 1 



+ 



^£sign(a-6)ir£-+i 



(6.5) 
(6.6) 



(6.7) 



where Z) = v^TlT~r is the diffusion constant. Through a similarity transformation S, 
we can bring this Hamiltonian to the form 



H' = SHS- 




q + q 



-1 



g + g 



-1 2 



aa 
+ 1 



a=0 



(6.1 



This similarity transformation does not have the simple local form of Eq. (|2.8|) . The 
argument goes as follows. We first notice that the Hamiltonian H' given by Eq. fl6.8|) is 
the exactly integrable U q SU(3) chain [36] and corresponds to the anisotropic version 
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of the spin 1 model introduced by Sutherland p7| , (see also Appendix C, where the 
connexion with the Hecke algebra is shown). On the other hand H' can be brought 
by an equivalent transformation to the four-parameter deformation of the SU(3) 
symmetric chain f38| 



H" 



SH'S 



-i 




10 

i+l 



+hoETET +l + f^E?Er +1 + hrE?E% x + f^E^E, 



t+i 



i \ - j-\aa jTiaa 

H 9 1^E % E i+1 



q-q 



-i 



5>ign(a - 



66 
i+l 



(6.9) 



We observe that if in Eq. (|6.9| ) we take /21 = /31 — fm — Q we recover Eq.(|6l 
The equivalence between the chain given by Eq.( |6.7| ) and the U q SU(3) chain 
is similar to the two-state model in the case of an asymmetric definition, where we 
got U q SU(2). The chain Eq. ( fi.tty is massive unless q — 1, which corresponds to the 
left-right symmetric case and the physical interpretation for these results is similar 
to the two species case of Section [|. The case where the process A + B — > B + A 
is suppressed is, as shown in Appendix C, again related to a quotient of a Hecke 
algebra and the eigenvalues of the Hamiltonian are the same (the degeneracies may 
be different) as for the U q SU(2) symmetric one given by Eq. ( |4.17| ). 

We shall take from now on q = 1 and parity conservation. We assume that the 
diffusion rates are 



^2,i(l,0) 

^2,l(2,l) 



^1,2(0,1) 
^1,2(1,2) 



W2,i(0,2) 
XD 



w 1}2 {2,0) = D 



(6.10) 



We consider now the annihilation processes 

with rate 101,2(0, 0) = iU2,i(0, 0) = olabE 
B with rate ^1,2(0, 2) = u>2,i(2, 0) = a^A-D 
A with rate 101,2(1, 0) = u>2,i(0, 1) = oibbD 

In this case, the Hamiltonian has a special form 

H = D(H + H 1 ) 




(6.11) 



(6.12) 



where 



L-l 



H = ^2 [^aaE^EI^ + 2a B BE? 2 Efl 1 



+ (A4 

+ R00 



E\lr + E 



77.11 77.22 
Ei E i+l 



77.22 77.11 
Ei E i+1 



22 

i+l 



77.IO 77.01 



77.01 77.10 
E t E i+1 



+ (E? + Ef) Ef* +1 

77.20 77.02 77. 
Ei E i+1 - E, 



02 77.20 
E i+ i 
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-A (if^r + E?E% X )] (6.13) 

H 1 = -J2[a AB (E^E^ + ETE^ 
i=i 

+a AA (jgf + £f i^) + ci_B_B (JSf£g a + EfE^)} (6.14) 

Using the arguments from Section one can check that the spectrum of the hermitian 
operator Hq coincides with the one of H/D. This Hamiltonian describes the pure 
reaction A + B — > 0, if we take 

A = a A A = a B B = (6.15) 
Let us find cases where the Hamiltonian is integrable. First we look at the 



U q SU(3) chain Eq. (|6.8|) with q = —1. If we add to the Hamitonian fields along 
the generators of the algebra, the system stays integrable since the supplementary 
terms commute with the integrable U q SU(3) Hamiltonian. We found only one case 
where this is possible. Taking 



A = 1 



OLAA — CHAB — GLBB 



we find 



L-l 

E 

i=i 



(6.16) 
(6.17) 



1 + E E ? E %i ~ E + 2 + 

The first three terms in Eq. ( |6.171 ) give a SU(3) symmetric Hamiltonian, correspond 
ing to q = — 1 in Eq. (|6.7|) . The last term with 



a 




(6.18) 



is an external field. This observation should allow us to study the critical properties 
of this chemical process. 

Another possibility occurs when A = 0. We can rewrite Ho in the following form, 
assuming o>aa = ®-bb 



L-l 

E 

i=l 



2a AA + OLAB 



2 ) + 



( Iqiaa - ®AB \ 



a i a i+l 



+ ^ + ^1 



T i T i+1 T i T i+1 Pi Pi+l Pi Pi+l 



where we have used the notations 
r+ = E 01 , r = E 10 



P 



E 



02 



P 



E 



20 



a 



E U ~E 
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This Hamiltonian commutes with E 2 and S° where 



i=l 



s° = E^° 

i=i 



(6.19) 
(6.20) 

(6.21) 
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For the special choice 



a A B = 2&AA 



(6.22) 



the Hamiltonian has the larger symmetry SU(2) <g>U(l) since H Q commutes also with 
£ which are defined by 



E 

i=l 



L 

i=i 



(6.23) 



-21 



using a + = E and a = E 

If we finally take olaa = 1, we have 



L-l 



£ [(<* 



a 



i+l 



Ti 7" 



i+l 



i+l Pi Pi+1 Pi Pi+1 



(6.24) 



i=l 



The chain Eq. ( |6.24 ) is integrable and massless as shown in Appendix C. It can also 
be related to t — J model [5*T] with J = 0. Moreover the spectrum of H (not the 
degeneracies) is that of a free fermion system. In Appendix C it is also shown that 
not only H but also H can be expressed in terms of a Hecke algebra which implies 
that the system is integrable. This remains valid even if we don't assume left-right 
symmetry. The quantum chain given by Eq.( 6.24j ) was discovered independently in 
Refs. ||53|1 and [Q. In the first reference this chain was found as a special case in a 
search of chains related to Hecke algebras. In the second reference it corresponds to a 
special case of integrable chains defined on nilpotent representations of the U q SU{2) 
quantum group with q = exp ^p. In the last reference the SU(2) eg) U(l) symmetry 
was observed. This triggered a further investigation by two of the authors to find 
the whole symmetry of the chain and obtain its eigenspectrum [55]. We would like 
to mention that for this system an exact result concerning the time-dependence of 



the total concentration < ca + cb > is already known | 3*9| . Moreover, if we keep 
the relation (|6.22j ) but let the parameter c^aa free, one can show that the spectrum 
(not the degeneracies) is the same as that of the XXZ model in a field, see Eq. ( |4.20[ ) 
with A = 1 — olaa and h — 1 — A, which corresponds to the Prokrovski-Talapov line 
of Fig. 1. Up to now we have found several choices for the constants aAB and o>aa 
for which Hq is integrable. We failed to find others, although two other classes of 
integrable models, namely the chiral Potts [HD| and the one with U q SU(2) , q 3 = 1 
symmetry where periodic or semiperiodic representations |4l|] are taken, have the 
required Z 3 symmetry. 

Much work is still to be done on the phase structure of the three-state models. 
One of the questions to be asked is the connection between the case Eq. ( |6.15|) , where 
®aa = 0, which is known to be massless, to the case ctAA = 1, which is massless as 
well. 

Finally, we mention that in Appendix C (see Eq. ( |CJ . 1 6| ) ) we give an example of a 
Z2 symmetric three-state chain (only A + A — > reaction) which is also related to a 
Hecke algebra. 
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7 Conditions on the reaction rates for the exis- 
tence of steady states 

Steady states are time-independent solutions of the master equation. In Section || we 



have seen that if the reaction rates satisfy the condition Eq. ( 2.12 ), then the empty 
lattice Eq. (gl3|) is a steady state. We now ask for the conditions on the rates such 
that the probability distribution 

k=l 

is a steady state, where f{/3) is a given arbitrary function. Our motivation is not 
only to understand better the structure of the ground states of the master equation 
and of the corresponding quantum Hamiltonians, but also to be able to address 
the question of the nonequilibrium behaviour of classical one- dimensional systems. 
This amounts to introduce a detailed balance condition for a prescribed equilibrium 
condition V, while in the previous sections we have always taken the empty lattice 
as the equilibrium configuration. Let us explain this point. 

Suppose that at temperature T we have a one-dimensional classical system defined 
by the probability function 

nm) = n #k - ^+0 (7.2) 

k=l 

If the /3's are site variables, the a's are link variables. 

In our study we are going to consider the two-state model (Ising model) where 

92(J3 k ) = 92(a k - a k+1 ) = exp {{-lf k /T) (7.3) 



and the three-state model (chiral Potts model [42]) where 

'2tt 



5"3(/3fc) = gz{ot k - afc+i) = ex P (cos (y/^ + / T ) ( 7 - 4 ) 
and is a fixed phase. We note the identities 

<fe(0)<72(l) = 1 , 03(O)03(1)<7 3 (2) = 1 (7.5) 



If the steady state is given by Eq. (|7.2| ), the master equation describes the evolution of 
the one-dimensional classical chain from P({j3};t = 0) = Po({/3}) to the equilibrium 
probability distribution P ({/?}; t = 00) = V({(3}). 

In order to find the conditions on the states, we first make a similarity transfor- 
mation Eq. ( |2.7|) , taking 



h {k) (a) = f 



1/2, 



a) 



W itm (a,/3) 



f(a + i)f((3 + m) 

\ JWW) ,m( } ( } 
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Then the condition to have Eq (|7.1| ) reads 



N-l 



w 0fi {a,P)= ^ wt, m (a,P) (7.7) 

faU f{a)f{(3) 

where the prime indicates that the case I = m = should be excluded. Taking into 
account the probability conservation condition Eq. (|2.4j) and making the change of 
functions 

wt, m (a, (3) = — — . , r 7.8 

f(a + l)f{f3 + m) 

we get the system of equations 

N-l 

Y! (7<,m(a, P) ~ le,mi a -t,P-m)) = (7.9) 

£,m=0 

We note that in Eq. ( |7.9| ) there is no longer any explicit reference to the function /. 
The whole /-dependence is in Eq. (|7.8| ). 

We now write down the independent conditions. To do so, it is convenient to 
introduce the following symmetric and antisymmetric combinations 

llmfaP) = Ji>,m(a,P)+JmAP' a ) 

llm{ a iP) = l/e,m(ai,P) -lm,l{Pi a ) ( 7 - 10 ) 



Now, for the two-state system the independent conditions coming from Eq. ( |7.9| ) are 
7i S ,o(0, 0) - 7l , (l, 0) = 7i,o(0, 1) - 7i S ,o(l, 1) = 7i,i(l, 1) - 7 i,i(0, 0) (7.11) 

7i%(l, 0) + 7i,o(l, 1) + 7o,i(l, 0) + 7o,i(0, 0) = 2^(0, 1) (7.12) 

Note that Eq. ( |7.11| ), which relates only symmetric combinations, does not, as ex- 
pected, contain the diffusion constant since it fixes the time scale only. The antisym- 
metric diffusion combinations appear however in Eq. (|7.12|) and thus play a dynamical 
role. We remind the reader that this property was already noticed in earlier sections 
when discussing Hamiltonians with quantum group symmetries (q ^ 1). 

For the three-state system with Z3 symmetry (which means that only the functions 
Wi t 2 and 102,1 appear) we have from Eq. (7.9) the following conditions 

72,i(0,0) = 72,(1,2) , 75,(1, 1) = 72,i(2,0) , 7 2,i(2,2)=72,i(0,l) (7.13) 



27^(1,0) = 7 2 a ,i(2,2)+7 2 a i(0,l) 
27^(0,2) = 7^(1, 1) +72,i(2,0) 

27^(2,1) = 7 2 a ,i(0,0)+7 2 a i(l,2) (7.14) 

We observe that the detailed balance equations do not determine all the rates. The 
remaining freedom can be used to fit the experimental data. We shall apply these 
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conditions on the rates in the next two sections. Before doing so, we briefly recall the 
questions of interest and some definitions relevant to critical dynamics. As already 
mentioned in the introduction, a classical Ising spin system has no intrinsic dynamics, 
because the Poisson brackets between the spin variables and the Hamiltonian vanish. 
In a real magnetic system, the spins are interacting with other degrees of freedom 
(phonons, impurities, . . . ) and those interactions are responsible for the dynamics of 
the spins. It is very difficult to implement fully the microscopic dynamics. Accord- 
ingly, one replaces the true dynamics by a fictious one given by a master equation 
for the probability that one spin configuration of the system is realized at a given 
time. The transition probabilities appearing in the master equation can in principle 
be computed from the full microscopic dynamics (via projective techniques). Un- 
fortunately, it is often not possible to perform this program, even for simple models 
f4*3| , Thus, the transition probabilities are chosen phenomenologically according 
to two criteria: i) to satisfy detailed balance, necessary to ensure that the desired 
equilibrium state is stationary, ii) to be qualitatively in agreement with the physics 
of the system. Once the dynamics is defined, the main problem is to find how the 
physical quantities will relax towards equilibrium. A quantity of particular interest 
is the order parameter (the magnetization for the Ising model). Typically, the order 
parameter relaxes very slowly in the vicinity of a second-order phase transition (criti- 
cal slowing-down). The dynamical scaling hypothesis assumes that the characteristic 
time r diverges as a power law of the spatial correlation length £; r ~ £ z , where z is 
called the dynamical exponent. 

The quantum chain formalism developed here aims to answer two types of ques- 
tions: 

1. Does the dynamical scaling hypothesis r ~ always hold ? 

2. What is the status of the "universality" for the dynamical exponent z ? 

For one- dimensional systems, the critical point is at zero temperature, T c = 0. 
For small T, the inverse spatial correlation lengths are 

= 2e~ 2 / T , two-state system 

r 1 = U^ + gj ^l) > three-state system (7.15) 
2 V#3(0) g 3 (0)J 

For the three-state system we shall only consider the ordinary Potts model, where 
= and the so-called superintegrable model |45| ) where <fi = 7r/6. This serves to 



illustrate how the dynamics depends on the equilibrium system. In the sequel, the 
deviation from the critical point will be parametrized in terms of the following masses 
fi ~ which are related to the temperature as follows 

fi = e~ 2 l T Ising model 

fi = e" 3/(2T) Potts model, = (7.16) 
fi = (r^l 1 ^) superintegrable model, <p = ti/6 
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for the three equilibrium models under consideration. In the next two sections we 
shall concentrate on the behaviour for small values of fi in the master equation and 
in the corresponding Hamiltonians. 



8 Critical dynamics for the Ising model 



We will write the Hamiltonian corresponding to the equilibrium distribution Eq. ( \l.3f j 
and the conditions Eqs. ( |7. 1 1| , |7.12| ) for the rates. We assume that the left-right 
antisymmetric combinations are zero such that Eq. (|7.12|) is trivially satisfied. To 
simplify the problem further, we choose the following solution of Eq. ( |7. 1 1| ) 

7i s ,o(0, 0) - 7i s ,o(l, 0) = 7i, (0, 1) - 7i,o(l, 1) = Ti,!^ l ) ~ 7i,i(0, 0) = (8.1) 
With this choice we get 
101,0(1, 0) = /iwi, (0, 0) , w 1>0 (l, 1) = ^i >0 (0, 1) , u; 1) i(l ) l) = /iVi(0,0) (8-2) 

This tell us that at T = 0(/x = 0) the master equation describing the critical dynamics 
of the Ising model reduces to the master equation considered in Sec. |] in which 
annihilation, coagulation and death processes were considered. The particular case 
with only annihilation was already noticed by Family aand Amar |32| . 
Taking take ti>i,i(l, 0) = 1 and by using probability conservation we get 



iw o> o(0,0) = /iV,i(0,0) + 2 Fl!0 (0,0) 
u> ,o(l,0) = l + wi, (0,0) +nw 1>0 (Q,l) 
tw 0)0 (l,l) = 2u> 1 ,o(0,l)+u> 1 , 1 (0,0) 

Using Eqs. ( |4.47| ), Q4.48D and (|8^) we get the Hamiltonian: 

L-l 



J.3) 



H 



-\ E [*f«SVi + + ^M+i + (i - a') (*? + 

i=l 

+ TOj i(0,0) («., • afrrf.i) 
+2^,0(0,0) (of E™! + 
+2 yU 1 /V,o(0,l) (ofE^+EFol 



T i+1 



•4) 



where 



A 
A' 



Wl>1 (0, 0) - (1 - /x)wi, (0, 1) + (1 - n)w 1)0 (0, 0) 



A - fiw 1>0 {0, 1) + /iV,i(0, 0) - (1 - 2//)wi )0 (0, 0) 



5.5) 



This result is very interesting and deserves a few comments. If ^1^(0,0) 7^ 0, even 
for /i = (when the equilibrium system is at the critical temperature T = 0) we 
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still have A' < A and consequently the time-dependent system is massive with an 
exponential fall-off in the correlation functions. This is to the best of our knowledge 
the first example of a system having an equilibrium second order phase transition 
but no critical slowing down. It is not obvious that this phenomenon generalizes into 
higher dimensions since only in one dimension, at T = T c the system is fully ordered. 

On the other hand if wi i o(0,0) = 0, but wi,o(0, 1) is non-zero, we would expect 
the relaxation time should be 

t" 1 -^ (8.6) 

The reason is simple and can be understood using perturbative arguments. We get 
terms of order 0{ji) in A', which couples to the U{1) scalars, of + of +1 , in the 
Hamiltonian. The terms present in the Hamiltonian which are of order 0(fx 1 ^ 2 ) occur 
in combinations like /i 1//2 wi j o(0, 1) which do not couple to U(l) scalars and should 
thus only contribute in second order. The contribution of order 0{n) can be only 
eliminated when a Z 2 symmetry is present in the problem. This implies that 

Wli0 (0,0) = Wli o(0,l)=0 (8.7) 

and by repeating the same argument we find 

r- 1 -/! 2 (8.8) 

However these perturbative arguments should be considered with some care since 
they apply to chains of a given length L. It might happen that in the large L limit, 
the leading term in the power expansions in \x has vanishing coefficients and we have 
to consider the next leading order term. To be more explicit, in a concrete calculation 
of the case Wi,o(0, 1) 7^ 0, iui,o(0, 0) = 101,1(0, 0) = for a given chain one obtains: 



.-1 



a L + b L fi + c L fi 2 + ... (8.9) 



In the large L limit ai must vanish (we end up in a massless system at fi — ) but 
the same can also happen to b^. It turns out that this is indeed happening since 



we can derive the exact result using the calculations of Ref . J3D| for the coagulation- 
decoagulation chemical process and find 

r 1 = £ (8.10) 

We would like to emphasize that the calculation of the coefficients aL,&z,,CL, ... can 
be difficult in general so that numerical estimates might be useful here. For the case 
^i,i(0, 0) 7^ and wi^O, 1) = wi^O, 0) the perturbative argument is correct and 



Eq. (|8.8|) stays valid as we can see from the exact result of Ref. |TL| . We conclude that 
in the Ising case the critical exponent z = 2 is universal. 
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9 Three-states critical dynamics 



Turning back to the three-states models, we now assume, in analogy to the two-state 
case, that the assymmetric rates in Eq. ( |7.14j ) are zero and we thus only have to solve 
Eq. (|7.13|) . We begin by considering the ordinary Potts model, with = 0. We find 



io 2 ,i(l,2) = /A/; 2 ,i(0»0) , w 2 ,i(l,l) = ^ 2 ,i(2,0) , w 2j i(2,2)=/itw 2> i(0,l) (9.1) 
This leads to the following Hamiltonian 



L-l 



,-.00 
i+l 



H = £[2/.VM(0,0)i^^ 

i=l 

+ (w 2A (l, 0) + ^ 2>1 (0, 1)) (E^Ell, + E?E™ X ) 
+ (w 2>1 (0, 2) + iiw 2>1 (2, 0)) (E^E 2 ^ + EfEl 
+ (w 2tl (2, 1) + ™ 2il (0, 0)) (E?E? +l + Ef E\\ x 
-^ 2il (0, 0) (El°E^ + E?Et* x + E?E™ X + E?Ef +x ) 
V'Vift 0) (E^Ell, + E?E% X + E?E$ X + E V2 E i 
V /2 ^,i(0, 1) {ETEll, + EfE 2 \ x + E?E™ X + EPS™, 

-^(1,0) (e™e°I x + e?eI 

-w 2>1 (2,l) (EfE^ + E^E] 
-w 2A (0,2) (E™E% 1 + E?°El 



-.10 
J i+1 



-.10 
i+l 

-.21 
J i+1 



,-.02 
i+l 



(9.2) 



In writing Eq. (|9.2|) , we have used the similarity transformation (|2.8| ) with h^(l)/h^(0) 
h^(2)/h^ k \0) = ^JJl. We are thus left, after choosing the time scale, with five free 
parameters. If we now take [i = in Eq. (|9.2j) , we get back to the chains studied 
in Section []. If we make the choice of coupling constants from Eqs. ( |6.10| , |6.12| ), at 
fi = we find that H coincides with the H of Eq. (|6.13|) . We can thus write 



H = H + H x 

where H x is in the case a AA = ®bb 

H x = Y[2^a AB ETE™ +1 + ^a AA ({l- p% P ' l+l + p^l 

i=l 

-^ 2 a AA ((e 10 + E 02 ) t EH, + El 2 (E 10 + E 02 ) t+i 
+ (£01 + E 20 )^ 2 ^ + Ef (E 01 + E 



(9.3) 



A+i 



(9.4) 
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and 

/ 1 \ 

p°= \=E 00 + E 22 (9.5) 

\o o l) 

Notice that the terms which couple to \iolaa keep the U(l) ® U(l) symmetry of the 
unperturbed problem. If this is the only perturbation present, we should thus get 

t- 1 ~ n (9.6) 



This result is similar to the behaviour observed for the Ising model (see Eq. ( |B.6|) ). 
We would like to stress that this result might be "naive" (see the discussion after 
Eq. (p|) in Sec. |). 

We now consider the superintegrable model, where = tt/6. This case is interest- 
ing because #3(2) = 1 and this will have important consequences. From Eq. Q7.13 ), 
we find 

w 2 ,i(l,2) = fi 3 w 2 ,i(0,0) , uj 2 ,i(M) =/i 3 ^ 2 ,i(2,0) , w 2>1 (2,2) = w 2>1 (0,l) (9-7) 
The Hamiltonian is 



H = £ [2^,1(0, 0)£^ 
i=i 

+ (w 2A (l, 0) + W2jl (0, 1)) (E^Ell, + E?E™ X ) 
+ (^(0, 2) + /^ 2>1 (2, 0)) (e™E% x + EfE™ x 



+ (w 2<1 (2, 1) + w 2A (0, 0)) (E^E^ + EfE\l x 
-^w 2A (0, 0) (El°E^ + + E?E% X + E?E$ X 

-^w 2A (2, 0) (E^E]l x + Efl x + E?El x + E?E™ X 



-w 2>1 (0, 1) (E™E% X + ETE^ X + E?E™ X + EfEf+ x 



10 

i+l 



-102,1(1,0) (EfEf^EfE 
-™ 2 ,i(2,l) (Sf^ 1 + ^ 2 ^ 



-IU2, 



,1(0,2)^^ + ^^)] (9 J 



Notice that for /i = we got a new term in Eq. Q9.8Q which did not exist in the 
Hamiltonians studied in Section ^| and which reads 



w 2 



,i(0, 1) (pfcrt+i + Pi + °i Pi+i + ViPtt-i) ( 9 - 9 ) 



where we have used the notation from Eq. ( |6.20| ). It is not known whether the 
Hamiltonian Eq. (|9.8|) for fi = is critical or not. Comparing the ordinary Potts 
and the superintegrable Potts cases we notice that the number of independent rates 
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is different. We also see, by comparing Eqs. (|9.2|) with (|9.9|) that ji 1 / 2 is replaced 
by /i 3//2 , consequently the dynamical exponents changes. We conclude that when 
moving from the Potts model (symmetry S3) to the superintegrable chiral Potts 
model (symmetry Z3), the exponents also changes. 

We have seen how the universal critical exponents can be predicted from con- 
sidering the structure of integrable Hamiltonians at /z = and appealling to per- 
turbation theory. Here again, we have implicitly assumed that all the six rates 
w 2,i(0, 0), W2,i(2, 0), u>2,i(0, 1), u>2,i(l, 0), u>2,i(2, 1) and w 2i i(0, 2) have finite, non-va- 
nishing limits as \x — > 0. Otherwise, one may create any value of some effective 



exponent like in Refs. 18, 43, 44 



10 Conclusions 

We have started our study asking ourselves what we thought is a simple question: 
is the present progress achieved in the understanding of one-dimensional integrable 
quantum chains useful for solving master equations describing the dynamics of classi- 
cal one- dimensional spin systems ? We find plenty of evidence for a positive answer. 
At same time, after finishing this long paper, we have the feeling that we are just at 
the beginning of a long path. 

Although initially we thought that our task would be just to use the available 
mathematical knowledge of integrable systems to find physical results, the physical 
problems brought a lot of feed-back into mathematics. Studying open chains with 
particular transitions rates lead us, to our surprise, automatically to new, hermitian 
and non-hermitian, representations of interesting associative algebras. We remind 
the reader that in equlibrium problems, quantum groups and associative algebras 
appear only through rather artificial boundary conditions. The physical applications 
of non-hermitian representations of associative algebras appear for the first time in 
our context. More about this subject will be published elsewhere. 

We would like to stress once more that although the phase diagram in the space 
of transition rates can be easily obtained from the Hamiltonian spectrum computed 
for example by using the Bethe ansatz, the calculation of nonequilibrium averages 
(which are not those normally occuring in equilibrium statistical physics) might pose 
formidable problems. 

For the physical understanding of react ion- diffusion processes or of simple critical 
dynamics we think that we went beyond the particular cases studied up to now in 
the literature. 

As the reader has certainly noticed while going through Sections |] and ^|, there is 
plenty of room for more work especially for three-states models, where we have only 
considered systems with Z3 symmetry and even in this case, the whole phase diagram 
is not yet completely known. 

Our experience with the study of properties of quantum chains for equlibrium 
purposes lead us to repeat analogous questions for nonequlibrium problems. For ex- 
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ample, in the simple model studied in Section |5| we have verified that finite-size scaling 
applies also to nonequlibrium situations. This opens the possibility of using standard 
numerical methods of matrix diagonalization of finite systems, for computing critical 
exponents. The solutions found in Section [7] for the detailed-balance equations should 
be useful for other dynamical processes than those discussed in Sections || and ||, in 
a better understanding of critical dynamics. 



Acknowledgements 

We would like to thank C. Gomez, V. Privman and G. Schiitz for useful discus- 
sions. We are most grateful to Klaus Krebs, Markus Pfannmiilller and Birgit We- 
hefritz for the excellent job they did critically reading the manuscript. F.A. thanks 
the Fundagao de Amparo a Pesquisa do Estado de Sao Paulo - FAPESP - Brasil and 
the Deutsche Forschungsgemeinschaft - DFG - Germany for support. M.D. and M.H. 
thank the Swiss National Science Foundation for support. V.R. would like to thank 
CERN, the University of Geneva, the Einstein Center of the Weizmann Institute and 
SISSA for the warm hospitality he has enjoyed and has enabled him to participate in 
this project. 

Appendix A. Eigenvectors of some special non-hermitian 
Hamiltonians 

We discuss the treatment of a certain class of non-hermitian Hamiltonians with the 
structure 

H = H + H 1 (A.l) 

where Ho is hermitian and has the same spectrum as H . Moreover we are interested 
in the case where the eigenspectrum is non-degenerated. Suppose the eigenvectors 
\ui) and the eigenvalues Aj of H 

H \m) = Xi \m) (A.2) 

are known and we want to find the eigenvectors of H. This problem is of interest for 
both analytical and numerical (finite-size scaling) calculations in reaction-diffusion 
processes (see Eqs. (4.11), (4.31) and (6.11)). In the basis Hi has only the 

following non- vanishing matrix elements 

H x \uj) = Gij ; i < j , i,j = l,...,N (A.3) 
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Thus the eigenvalues of H are again Ai, . . . , \n- Let \vi) be the eigenvector of H 
corresponding to the eigenvalue Aj. It can be written as 



(A.4) 



with An = 1 and the graphical rule to compute the Ay(i > j) is obvious if we give 
the first few 



A21 
An 



G 



12 



A2 — Ai 
A4 — A3 



^32 
^42 



G 



23 



A3 — A2 



A 



G 



13 



G10G 



12^23 



31 



+ 



A3 — Ai (A3 — A2XA3 — Ai) 
C23C34 



A4 ~ A2 (A4 — As)(A4 — A2) 



+ 



G13G 



13*^34 



+ 



GtoG 



12^24 



A 4 - Ax (A 4 - A 3 )(A 4 - Ax) (A 4 - A 2 )(A 4 - Ai) 
G12G23C34 
(A4-A 3 )(A 4 -A 2 )(A4-A 1 ) 



(A.5) 



Appendix B. The seven- vert ex model 



As an example, we shall consider in more detail one of the models related to a Hecke 
algebra (see Appendix C) in order to obtain the Boltzmann weights of an associated 
two-dimensional vertex model and thus prove its integrability. 

The model we consider is the one introduced in Section ^ which describes diffusion 
and pairwise annihilation of A particles, see Eq. (|4.8| ). We choose the special tuning 
which makes A = and take the functions ip( k \a,/3) = 1 where ot,/3 = 0,1. The 
Hamiltonian is given by 

H = D(H + H 1 ) (B.l) 



where 



Hn 



L-l 

E 

i=l 
L-l 



qE?E™ x + -EfE^ + 1 (at - 1) + f <i - 1 



2q 



01 Z7.01 
t+1 



and 



n 



-n £ e^e\ 
1=1 

wi,i(0,0) 



9 



\ 



IB.2) 



(B.3) 



Doing the canonical transformation i£. 



— 1) -Ej , only at even sites i, this 



Hamiltonian in the a z basis takes the simple form 

H 



L-l 

£>E e * 

i=i 



(B.4) 



41 



where 



li <g> ... <g> lj_i <g> 



/ 







1 

'/ 
1 






Q 
Q 




n 




1 +L J 



H+2 



<8> . .. 



(B.5) 



and i = 1,2, . . . ,L — 1 and lj are 2x2 unit matrices attached to the site z. 

The Hamiltonian Eq. ( |B.4|) is known to be integrable through a Jordan- Wigner 
transformation (see Section we think however that the approach given here gives 
not only new insight into the problem but is of a larger validity. 

We can show that the above matrices ei,i = 1,...,L satisfy the Hecke algebra, 
for arbitrary values of Q 







L e i> e j\ 



'- ej±iejej±i — ei±\ 
>2 



(B.6) 



q + q 



This is the first example we know of where nonhermitian (take q real) representations 
of the Hecke algebra appear in physical applications. The hermitian case Q = 
corresponds to the quantum chain introduced by Saleur [|6 . 

Due to the algebraic properties Eq. (|B.6| ) we can construct an associated two- 
dimensional vertex model having a row-to-row transfer matrix depending on the 
spectral parameter 9. There transfer matrices will satisfy the Yang-Baxter equations 
|47| which implies that they commute among themselves for different values of the 
spectral parameter. 

In order to obtain the configuration and the Boltzmann weights associated to this 
vertex model, we need the spectral parameter dependent matrix Ri(9),i = 1,2, . . .. 
This is found by the Baxterization procedure [BTJ for Hecke algebras, namely 



sinh 6 sinh(?7 

-e, ; + 



which gives us 



Ri(0) = h 



Li-IA 



sinh 7/ 



sinh 7] 



( sinh (77 



V 



sinh?] 



q 



(B.7) 



e sinh 77 e v sinh 9 
e~ v sinh 9 e e sinh r] 



Q sinh 9 \ 



sinh (r/ + 9) ) 



(B.8) 

The relations Eq. ( |B.6| ) imply that these matrices satisfy the spectral parameter 
dependent braid group relations 



R i (6)R i ±i(6 + 6')R i (9' 
'R t (9),R J (9') 



Ri±i(9')Ri(9 + 9')Ri±i(9) 
, \i-j\>2 



(B.9) 
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which are equivalent to the Yang-Baxter equations. 

The Boltzmann weights Sf^ of the vertex configuration labelled by (k, £, m, n) in 
the associated vertex model can be obtained from the relation 

RiW) = Sf/mll ® • ■ ■ ® 1<-1 ® ^ ® ^ ® li+2 ® • • • (B.10) 

This implies that the vertex model associated to Eq. ( B.l ) is a seven vertex model. 



If we denote by an index zero a down (or left) arrow and by an index one an up 
(or right) arrow, the vertex configurations with their Boltzmann weights are given in 
Figure 2. We also show in this figure the corresponding chemical processes related 
to each vertex. The vertices 1 to 4 correspond to no reaction, vertices 5 and 6 to 
diffusion to the right and the left and vertex 7 to the pair annihilation process. The 
derivative of the logarithm of the row-to-row transfer matrix, with these Boltzmann 
weights and evaluated at 9 = 0, gives back the Hamiltonian Eq. (p.l|) . 



Appendix C. Hecke algebra and reaction-diffusion 
processes 

In this appendix we define the Hecke algebra and give some examples of Hamiltonians, 
related to dynamical processes, which are representations of this algebra. 

The Hecke algebra is an associative algebra with generators ej (i — 1, . . . , L — 1) 
satisfying the relations 

eje i±1 ei - e» = e i ± x e i e i ± x - e i±1 (C.l) 
[e h ej] =0 ; \i-j\ >2 (C.2) 

(q + T 1 ) e* (C3) 



el 



where q is a complex parameter. 

One of the main features of the above algebra is related to the fact that a spectral- 



dependent R{u) matrix, satisfying the Yang-Baxter relations |47|], can be constructed 
in a standard form. As a consequence of this procedure, also called "Baxterization" 
2jJ , the Hamiltonian J2i=i e i> as weu as &s associated vertex model, has an infinite 



number of conservation laws and we expect, in general, its exact integrability. In 
Appendix B we give an example of this "Baxterization" procedure and derive the 
associated vertex model for one of the chains considered in this paper. 

It is important to stress here that distinct Hamiltonians satisfying the same Hecke 
algebra may correspond to different representations of the algebra. They will share, 
apart from degeneracies (which may be zero) the same eigenenergies. These chains 
will have a massive or massless behaviour depending on the value of q. In particular, 
for q real they will always have a massive behaviour except for q — 1, where they will 
be massless. 

The ej's appearing in quantum chains quantum chains satisfy beyond Eqs. ( |C.1| - 
]3|) additional relations which define quotients of the Hecke algebra. Obviously, 
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quantum chains obeying the same quotient have more coincidences in their spectra, 
see |23] for details. A well-known quotient is the Temperley-Lieb algebra |^8[ defined 
by QOJ) and 



1,2, 



(C.4) 



A realization of this algebra is given by the XXZ chain with surface magnetic fields 
(invariant under the U q SU (2/0) quantum group) and the quantum Potts chain with 
free ends |48|, |49fl . In Sections [| and |^ we have used the notation U q SU{n) for 
U q SU(n/0). With the new notation U q SU (n/m) with m/0 corresponds to quantum 
superalgebras p0| .) Another less known quotient is defined by ( |C.1| - |U.3| ) and the 



additional relation 



eie i+ 2)e i+1 [q + q - eA [q + q 



e;+2 







1,2, 



(C.5) 



invariant under the 



which has a realization in the two-colors Perk-Schultz model p2 
quantum group U q SU( 1/ 1 ) . 

In this paper we show that several quantum chains, related to chemical processes, 
are realizations of the Hecke algebra (C.1)-(C.3). As a general outcome from our 
analysis we verifyed that as long as only diffusion and interchange of particle processes 
(the number of particles in each species is conserved separately) are allowed, these 
realizations arise quite naturally In the cases where other processes are also allowed 
the chain in general will not satisfy the Hecke algebra. However, as we shall see, 
for certain processes and special tunings of transition rates, the corresponding chains 
will turn out to satisfy the Hecke algebra again. 

Although we can generalize our results to an arbitrary number of different chemical 
species, following the line of this paper we will consider here only the cases where we 
have, beyond vacancies, 1 species (A) or 2 species (A and B). 

Let us consider initially the cases where we have only particles and vacancies, see 
Section f|.l. If we only allow the diffusion process with the transition rates satisfying 
Eqs. ( |4.9D , ( [4. lip we obtain the XXZ chain with anisotropy A = (q + q~ 1 )/2 



H/D = H 



L 

i=i 



(C.6) 



where 



&i + 0; 07-4-1 + - I P + ~ I Gi &i 



"i u i+l 





(C.7) 



We can verify that these matrices are just a 2 L -dimensional representation of 
the Hecke algebra. Moreover, this Hamiltonian satisfies the quotient defining the 
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Temperley-Lieb algebra flG.4j) . It is important to stress that the terms proportional to 
erf appear naturally in ( |C.7| ) due to the diffusion mechanisms. These terms, although 
not present in a periodic chain, are crucial in order to generate the Temperley-Lieb 
algebra. The particular chain ( |U.7|) is invariant under the quantum group U q SU (2/0). 

Another example, also with only particles and vacancies, appears when beyond 
diffusion we also allow the annihilation process A + A — > + (with rate wi t i(0, 0)). 
In this case, by making the special tuning 



wi,i(0,0) = w 1) i(0,l)+w 1 , 1 (l,0) 
and choosing the functions ip( k ' (£, m) = 1 we obtain 



H/D 



L-1 

1=1 



(C.9) 



where 



i+1 



^ Z?10 77.01 




+~(crf-l) + - (a. 



2q 



\ u i+i 



L-1 



and 



wi,i(l,0) 



D 



^i, 1 (o,iKi(i,o) 



(C.IO) 



(c.n; 



H ^ 1^,1(0,1) 

These matrices satisfy the Hecke relations (C.1)-(C.3) and are also generators of the 
quotient defined by Eq. ( |C.5| ). This is the first example we are aware of in which non- 
hermitian (q real) representations of the Hecke algebra appear in a physical context. 
In Appendix B we derive the vertex model associated to this chain. 

Let us now consider some cases where we have two types of particles (A and 
B). If we consider only the process of diffusion and interchange of particles, as we 
saw in Section || by choosing the diffusion rates Eq. (|6.4j ) we obtain the anisotropic 
Sutherland model |!6| 

H/D=J2e l (C.12) 



where 




a,b=0;ajtb / \ 1 / a=0 



-, £ sign(a - b)E?E* r 

1/ a,b=0;a^=b 



(C.13) 



which again satisfy the Hecke algebra. The above Hamiltonian is invariant under 
the quantum group U q SU(3/0) (we do not give the corresponding quotient here). 
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Let us now consider the case where beyond the above processes we also include the 
annihilation A + A —>■ + (rate Wi,i(0, 0)). If we use the relation (6.5) between the 
diffusion rates and the condition 



we obtain 



where 



wi,i(o,o) = r x + r fl 



H/D' = J2 e % 



Ir 01 R 10 + n E w R m +-E 02 E 20 



(C.I4) 
(C.15) 
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+qE?E% 1 + -E?E& X + qEfEfa + QE^E^ 
+Q (E^Ell, + E™Ef +l + E^EJl,) + \ (E?E™ X + E?E™ +1 + E?E„ 
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i+l 



9, 



and 



T R /T L , D' = JY R T L 



(C.I6) 



(C.17) 



The (nonhermitian!) operators e; also satisfy the Hecke algebra. In fact we verified 
that this is a property of ( |C.16|) - (|C.17|) for arbitrary values of Q. For Q = the 
Hamiltonian reduces to the three-color Perk-Schultz model [^2] (which has the quan- 
tum superalgebra U q SU (2/1) as symmetry) and it is also related to a special point 
of the t — J model pnj where exact integrability takes place [jj^ • 

Let us return to the case where we do not have annihilation, see Eqs. (|C.9| , |C.lTD . If 
we now forbid the process where the particles interchange positions (A+B <-> B + A) 



wi,2) = W2 i(2,r 



we obtain 



with 



L-l 



H/D' = £ e, 



(C.18) 
(C.19) 
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E K<i + EfE^r) + qEi™ (EH, + 



.0=1 



--[Ej' + Ef 
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i+l 



1 L 



(C.20) 
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This model was introduced in |53| and the generators satisfy the Hecke algebra. 
Moreover we verified that like the U q SU (2/0) model ( p.6|JC7| ) the matrices e; are 
also the generators of a 3 L -dimensional Temperley-Lieb algebra. The symmetry of 
the chain ( |C.20| ) is known and described in Ref . . 

If we now include (see Eq. (|6.11|) ) the annihilation processes A + B — > + 0, 
A + A — > $ + B, B + B — > + A with the reaction rates related in the special way 

w 2>1 (l, 0) = w 1>2 (2, 0) = w 2>1 (2, 0) = ^,2(1, 0) = T L (C.21) 

Wl 2 (0, 1) = w 2)1 (0, 2) = w 1>2 {0, 2) = tu 2>1 (0, 1) = T R (C.22) 
w 1>2 (0, 0) = w 2A (0, 0) = T R + T L (C.23) 



we obtain 



where 



q = \JTl/T r , D' = ^T L T R (C.24) 
H/D' = H /D' + Hx/D' = J2 e, , = e° + e] (C.25) 



i=l 



a=l\ Q Q / 

+ U+-) ; < = 1,...,L-1 (C.26) 



and 



e. 



(-E^Ell, + qE^E^ + qE^ + UfE] 



02 zt.12 
+ 1 



+ [q + I j + (? + ^) (C27) 

The Hamiltonian H /D' = Y.i=i e i was introduced in (see also Ref. for the 
case g = 1) as a 3 L -dimensional representation of the Hecke algebra. We verified, 
guided by the physical processes of diffusion and annihilation, the the operators e« 
(i = 1, . . . , L — 1) also satisfy the same algebra. Morevover, we checked that both 
e° and e; satisfy the quotient (|C.5| ) like the U q SU (1/1) chains, for the underlying 
quantum symmetry of ( |C.25| ) see Ref. |)55| . 

Figure captions 



Figure 1: Phase diagram of the XXZ quantum chain with anisotropy A in a magnetic 
field h. I is the massive (frozen) ferromagnetic phase, II the massless phase which 
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is commensurate for h = and incommensurate for h 7^ and III is the massive 
(frozen) antiferromagntic phase. The line A = 1 — h corresponds to a Pokrovsky- 
Talapov transition. 

Figure 2: Boltzmann weights of the two-dimensional vertex model associated to the 
Hamiltonian Eq. ( |B.1|) . The first column gives the weights, the second one the arrow 
diagram and the third one the relationship between the vertices and the elementary 
reaction-diffusion processes. 
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